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Abstract 


The  response  of  the  neutral  atmosphere  to  momentum  forcing  at  F-region  heights 
(about  400  km)  is  investigated.  The  appropriate  form  of  the  momentum  equation  as  it 
pertains  to  the  observed  local  acceleration  of  the  wind  at  a  specified  location  is  developed. 

Experiments  were  performed  using  data  generated  by  a. National  Center  for  Atmospheric 

—  — . —  — 

Research  (NCAR)  thermospheric  general  circulation  model  (TGCM)  The  acceleration  of 
the  neutral  wind  at  several  locations  and  times  were  compared  with  the  sum  of  the  forcing 
terms  of  the  momentum  equation.  The  appropriate  forms  of  the  equation  for  stationary  as 
well  as  moving  observers  were  used,  the  latter  being  applied  to  satellites  orbiting  at 
thermospheric  altitudes.  The  momentum  equation  was  analyzed  to  determine  which  forces 
could  be  calculated  solely  from  the  data  typically  gathered  by  the  Dynamics  Explorer  (DE- 
2)  satellite  Both  synthesized  and  actual  satellite  data  were  used  in  the  various  experiments. 

Techniques  and  software  were  developed  to  assist  in  the  data  handling  and 
subsequent  analysis.  The  results  ranged  from  excellent  agreement  between  the  observed 
acceleration  and  total  forcing  for  the  case  of  a  fixed  observer  when  archived  TGCM  data 
was  used,  to  virtually  no  correlation  between  measured  momentum  forcing  and  observed 
acceleration  in  the  case  of  the  satellite- borne  measurements.  It  was  found  that  the  best 
results  were  obtained  when  the  effect  of  ion  drag  on  the  neutral  wind  was  at  a  minimum. 
This  occurred  during  the  evening  hours  (when  ion  densities  were  lowest)  and  at  low 
latitudes  (away  from  polar  electric  fields).  A  detailed  analysis  on  measurement  uncertainty 
and  its  potential  effect  on  the  results  is  also  presented. 
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Chapter  1  Introduction 


1.1  Statement  of  Problem 

The  objective  of  this  thesis  is  to  study  the  response  of  the  neutral  thermosphere  to 
momentum  forcing  inputs  and  verfy  that  the  measured  response  is  consistent  with  the 
momentum  equation. 

To  accomplish  this  object  I  have  made  extensive  use  of  the  output  of  a 
thermospheric  general  circulation  model  (TGCM)  to  simulate  the  evolution  of  the  neutral 
winds  at  thermospheric  heights.  This  data  has  been  examined  extensively,  and  detailed 
comparisons  were  made  between  the  forces  acting  at  specific  locations  and  the  observed 
accelerations  at  those  locations. 

1.2  Thesis  Format 

This  thesis  consists  of  five  chapters  and  four  appendices.  Chapter  2  is  a  brief 
review  of  the  momentum  equation  as  it  pertains  to  the  thermosphere.  The  lagrangian  and 
eularian  forms  of  the  momentum  equation  are  presented  as  is  the  application  of  the  latter 
form  to  the  present  problem.  The  ion  drag  and  viscous  drag  terms  are  described  as  being 
the  respective  major  and  minor  driving  forces  in  this  atmospheric  region. 

Chapter  3  describes  the  major  tools  that  were  used  for  this  study.  The  TGCM  and 
VSH  models  are  briefly  described  and  several  examples  of  their  output  are  presented.  The 
working  data  base  from  which  the  majority  of  the  computations  was  made  is  also  described 
here.  Finally,  a  description  of  the  Dynamics  Explorer  (DE)  on  line  data  base  at  the 
University  of  Michigan  is  presented  along  with  some  information  as  to  how  it  was 
specifically  applied. 

The  main  core  of  the  work  is  presented  in  chapter  4.  It  is  in  this  chapter  where  the 
momentum  equation  is  broken  down  into  its  component  parts  and  detailed  comparisons 
between  the  forcing  and  acceleration  are  made.  The  process  begins  by  looking  at  archived 
TGCM  output  and  comparing  the  acceleration  of  the  neutral  wind  to  the  sum  total  of  all  the 
forces  acting  at  one  specific  location  over  a  24  hour  period.  Next,  the  same  analysis  is 
performed  for  this  and  other  locations  using  the  data  from  the  working  data  base. 
Following  the  fixed  location  study  the  same  sort  of  analysis  was  performed  for  the  case  of  - 
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a  simulated  polar  orbiting  satellite.  The  chapter  ends  with  a  study  of  thermospheric  forcing 
based  on  observerations  made  by  a  polar  orbiting  satellite  using  actual  data  from  the  DE 
data  base.  Results  and  conclusions  are  presented  in  chapter  5. 

The  appendices  contain  useful  information  that  was  not  necessarily  appropriate  for 
inclusion  in  the  main  body.  Appendix  A  is  a  detailed  look  at  the  concept  of  the  total 
derivative  and  its  applicability  to  the  present  problem.  As  a  thorough  understanding  of  this 
concept  is  crucial  to  this  study,  I  would  strongly  recommend  reading  Appendix  A  before 
beginning  chapter  2.  Appendix  B  contains  the  techniques  and  formulas  used  to  compute 
the  uncertainties  in  the  various  computations.  Appendices  C  and  D  contain  tables  and 
graphs  of  the  simulated  and  actual  satellite  data  that,  though  important  for  computational 
purposes,  was  not  crucial  for  the  discussion  in  the  main  body  of  the  thesis. 
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Chapter  2  The  Momentum  Equation 


The  evolution  of  the  motion  of  any  portion  of  the  atmosphere  is  ultimately  governed 
by  the  fundamental  equations  of  motion.  For  any  individual  control  volume  (or  parcel)  of 
air  the  motion  can  be  described  by  the  following  equation  as  presented  by  Holton  [1979]: 


dU 

dt 


=  -2Q  x  U  -  -jjVP  +  g  +  Fr 


(2-1) 


In  the  above  expression  the  term  on  the  left  hand  side  is  the  acceleration  of  the  control 
volume  following  the  motion  of  the  fluid  (the  airstream  in  this  case).  The  first  term  on  the 
right  hand  side  is  the  Coriolis  force  that  arises  because  of  the  accelerated  nature  of  the 
earth’s  rotating  reference  frame.  The  second  term  on  the  right  is  the  pressure  gradient  force 
common  to  all  fluid  motions.  The  third  term  is  the  acceleration  due  to  gravity  and  always 
acts  downward  (toward  the  origin  in  an  earth-centered  coordinate  system).  The  last  term 
on  the  right  hand  side  (Fr)  is  the  sum  of  forces  due  to  friction  with  adjacent  flows.  It  is 
expressed  in  generic  form  here  because  its  specific  form  is  unique  to  the  portion  of  the 
atmosphere  under  consideration.  The  form  of  this  term  appropriate  for  this  study  will  be 
discussed  shortly.  During  this  study  we  consider  the  two-dimensional  form  of  the 
momentum  equation  (2-1)  for  neutral  species.  Thus,  the  eastward  (u)  and  northward  (v) 
horizontal  components  of  the  neutral  wind  will  alone  be  considered;  so  the  g  term  in  (2-1) 
does  not  come  into  play. 

The  coordinate  system  used  in  this  study  is  the  one  most  commonly  used  to 
describe  locations  and  motions  in  the  geocentered  system;  namely,  the  spherical  coordinate 
system.  The  form  of  (2-1)  appropriate  for  a  rotating  spherical  coordinate  system  is 
presented  by  Holton  [1979]as  follows: 


dU  _  f  du_  _  uv  tan  0  \  f  dv  »2tan  0\ 
dt  V  dt  r  7  V  dt  r  J3 


(2-2) 


where:  0  =  latitude 

X  =  longitude  (below) 
r  =  radius  of  earth  plus  altitude 
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du/dt  =  instantaneous  zonal  acceleration 
dv/dt  =  instantaneous  meridional  acceleration 


The  curvature  terms  in  (2-2)  arise  because  of  the  positional  dependence  of  the  unit  vectors  i 
and  j.  The  individual  component  momentum  equations  then  become: 

1 


du 

dt 


prcos  0  dX 


dP  .  c  .  r  ,  uv  tan  0 
+  fv  +  F.  +  - 7 - 

*  k 


(2-3a) 


dv  _  1  3P  ,  ^ 

dt  "  pr  30  fu  Fy  r 


(2-3b) 


2.1  Thermospheric  Form  of  the  Momentum  Equation 

The  form  of  the  momentum  equation  for  the  region  of  interest  (ionospheric  F- 
Layer)  is  found  by  substituting  the  appropriate  frictional  forces  into  (2-2).  These  frictional 
forces  include  the  viscous  drag  force  md  the  ion  drag  force. 


2.1.1  Viscous  Drag  Force 

The  viscous  force  is  due  to  the  action  of  vertical  diffusion  in  the  presence  of  vertical 
velocity  shears.  The  expression  describing  the  magnitude  of  the  viscous  drag  force  is  also 
described  by  Holton: 


'I  32u 
x  P  9z2 


(2-4a) 


F 


y 


d2v 

P  3z2 


(2-4b) 


where:  q  =  coefficient  of  dynamic  viscosity 

p  =  mass  density  of  neutral  atmosphere 

In  the  above  expression  the  q  term  is  a  slowly  varying  function  of  temperature,  while  the 
mass  density  varies  significantly  with  altitude.  The  ratio  q/p  is  often  called  the  kinematic 
viscosity,  and  also  has  a  marked  altitude  dependence.  Thus,  we  would  expect  the  viscous 
drag  to  work  efficiently  at  the  F-region  altitudes.  Ironically,  this  efficiency  quickly 
removes  the  second  derivative  in  the  velocity  profiles  and  subsequently  reduces  the 
magnitude  of  the  viscous  drag  force.  As  a  result  the  viscous  drag  force  has  a  relatively 
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minor  role  in  the  horizontal  momentum  equation.  The  next  chapter  quantitatively  compares 
the  viscous  drag  with  the  other  terms. 

2.1.2  ’on  Drag  Force 

The  F-region  of  the  ionosphere  contains  a  high  density  of  positively  charged  ions 
mostly  in  the  form  of  singly  ionized  atomic  oxygen.  Because  of  the  net  charge  on  the  ions 
they  can  be  accelerated  in  the  presence  of  electric  fields.  This  acceleration  is  completely 
independent  of  the  i..:>tion  of  the  surrounding  neutral  gas,  and  a  velocity  shear  between  the 
two  constituents  is  almost  always  present.  As  ions  collide  with  the  neutral  particles  a 
frictional  drag  force  is  imparted  to  the  neutrals  according  to  the  following: 

n . 

As  can  be  seen  from  (2-5)  the  magnitude  of  the  ion  drag  force  depends  not  only  on  the 
velocity  shear,  but  on  the  relative  density  between  the  ions  and  neutrals  as  well.  This  fact 
creates  a  local  time  dependence  in  the  ion  drag  as  the  ion  population  goes  through  a 
pronounced  diurnal  variation.  The  momentum  transfer  collision  frequency  can  be 
computed,  and  the  techniques  for  doing  so  are  included  in  Appendix  B.  As  we  will  see  in 
the  next  section,  the  ion  drag  force  is  a  major  player  in  determining  the  neutral  gas 
circulation. 


(2-5a) 

(2-5b) 


2.2  Momentum  Equation  in  Eulerian  Form 

Upon  substituting  the  viscous  and  ion  drag  forces  into  (2-3)  the  momentum 
equations  following  the  motion  of  the  air  (the  so-called  lagrangian  form)  become: 


du 

dt 


1  ap  +  fv  +  n 


uv  tan  <j) 


prcos  0 


F^T+Vui-un)  +  — r 


t2-6a) 


dv  _ _ 1  dP  _  f  *1  d2  v 

dt  pra<j>  p  az2 


+  V  (v.  -  v_)  - 
m'  i  n/ 


u2tan  <f> 


(2 -6b) 


An  alternative  to  the  lagrangian  form  above  is  the  eulerian  representation,  which  relates 
accelerations  following  the  motion  to  those  observed  at  a  fixed  location  as  follows: 


dQ 

dt 


=  ^+V.VQ 


(2-7) 
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Using  this  expression  for  the  total  derivative  we  can  express  the  momentum  equations  in 
the  following  eulerian  form: 


9u 

9t 


1 


9P  ......  *1  9  u  ,  -  ,  .  uvtan<)) 


prcos  <p  dX 


+  tw*p-^2+VJU,  -“«)  + 


-  V  Vu 


9v 

9t 


J_9P  f  .  H.  92 v  _  ,  ,  u2tan0 

pr  9$  fu+  P  9Z2  +  Vn/vi  v n  r 


-  V  •  Vv 


(2-8a) 

(2-8b) 


We  will  make  extensive  use  of  this  form  of  the  momentum  equation  both  in  the  case  of  a 
fixed  observer  and  in  the  case  of  an  observational  platform  that  is  moving  with  respect  to 
the  established  coordinate  system  (ie  an  orbiting  satellite).  A  complete  development  of  the 
physical  meaning  of  the  total  derivative  and  how  it  relates  to  stationary  as  well  as  moving 
observers  is  presented  in  Appendix  A. 

As  it  stands  the  momentum  equation  gives  a  complete  description  of  the  forces 
acting  at  a  specified  location,  but  it  gives  no  information  as  to  which  of  the  forcing  terms 
are  the  most  significant.  Figure  2.1  is  a  stacked  breakdown  of  the  individual  forcing  terms 
during  a  24  hour  period  at  location  (  47.5,  0  ).  From  this  figure  one  can  see  that  the 
pressure  gradient  and  the  ion  drag  are  clearly  the  most  important  contributors  to  the  local 
forcing.  The  forces  depicted  on  the  diagram  were  computed  using  a  thermospheric  general 
circulation  model  described  in  the  next  chapter. 
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Chapter  3  Tools  used  for  this  Study 


This  study  was  largely  computer  based,  both  in  the  area  of  modeling  and  in  data 
retrieval  and  analysis.  The  modeling  was  accomplished  primarily  with  the  National  Center 
for  Atmospheric  Research  (NCAR)  thermospheric  general  circulation  model  (TGCM).  The 
TGCM  and  its  diagnostics  processor  was  used  to  build  a  large  working  data  base  from 
which  momentum  forcing  calculations  could  be  made.  In  addition  to  the  TGCM  a  vector 
spherical  harmonic  (VSH)  decomposition  model  was  used  to  supplement  TGCM  data. 
Finally,  a  large  data  base  of  actual  measurements  from  the  Dynamics  Explorer  2  (DE-2) 
was  available  and  used  for  comparative  study.  This  chapter  briefly  summarizes  the  models 
and  the  content  of  the  working  data  base. 

3.1  Thermospheric  General  Circulation  Model  (TGCM) 

The  momentum  forcing  aspects  of  the  TGCM  as  well  as  a  description  of  the 
diagnostics  processor.is  discussed  by  Killeen  and  Roble  [1984],  This  model 
si  nultaneously  solves  the  continuity,  momentum,  and  energy  equations  to  compute  the 
wind,  temperature,  and  composition  for  all  model  grid  points  at  each  hour  of  universal 
time.  The  grid  is  defined  horizontally  by  the  5  degree  latitude  and  longitude  points  and 
vertically  by  24  constant-pressure  surfaces  with  1/2  scale  height  separation.  The  model  can 
accept  several  geophysical  parameters  as  input  and  is  thus  very  useful  for  studying  the 
evolution  of  the  thermosphere  under  a  variety  of  conditions.  The  model  is  typically  run 
until  a  steady  state  is  reached  and  a  diumally  reproducible  history  file  is  created.  The 
history  file  used  in  this  study  was  created  using  the  following  geophysical  parameters  as 
inputs: 


Hemispheric  power: 

11  GW 

Cross  Cap  Potential: 

60  kV 

F10.7: 

243 

Equivalent  Day: 

79355 
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The  diagnostics  processor  was  developed  to  provide  a  means  for  making  detailed 
analysis  of  the  individual  forcing  terms  of  momentum  equation  for  each  location  and  at  each 
universal  time.  The  processor  reads  the  TGCM  history  file  as  its  primary  input,  but  also 
accepts  some  of  the  same  inputs  that  were  used  in  the  model  run.  Earlier  TGCM  versions 
[Dickinson,  et  ai,  1981]  include  (among  other  things)  the  ion  drag  tensors  used  to  compute 
the  ion  drag  force.  Later  versions  of  the  model  (including  the  one  used  here)  computed  the 
ion  drag  tensors  within  the  code  itself.  The  processor  was  run  on  the  history  file  described 
above  to  create  the  working  data  base. 

3.2  The  Working  Data  Base 

For  each  horizontal  grid  point  the  working  data  base  initially  included  the  neutral 
and  ion  winds,  neutral  and  ion  temperatures,  neutral  and  electron  densities,  height  of  the 
pressure  level,  ion  drag  tensors,  and  the  mass  mixing  ratios  of  molecular  and  atomic 
oxygen.  Later  the  individual  forcing  terms  were  also  added,  which  includes  the  ion  drag, 
pressure  gradient,  Coriolis,  viscous  drag,  the  momentum  advection,  and  the  sum  total  of 
the  individuals.  Only  one  pressure  level  (Z=2)  was  considered  for  this  study,  the  pressure 
levels  being  defined  by  the  following  simple  expression: 

Z  =  logeOVP)  (3-1) 

where  Po  is  equal  to  the  reference  pressure  of  50  jiPa.  From  this  expression  we  see  that 
integral  values  of  Z  are  separated  by  a  factor  or  1/e  (ie  one  scale  height).  The  altitude  of  the 
Z=2  pressure  level  was  approximately  375  km  for  this  study.  Even  though  only  one 
pressure  level  was  considered,  the  volume  of  the  retrieved  data  all  but  filled  all  allocated 
disk  storage  space.  To  make  a  more  manageable  working  data  base,  only  the  grid  points 
corresponding  to  latitudes  greater  than  or  equal  37.5  degrees  (north)  were  retained  for  this 
study.  The  intact  data  base  was  transferred  to  magnetic  tape  for  future  availability. 

Each  variable  in  the  data  base  was  contained  in  its  own  data  file  and  each  had 
identical  file  structure.  The  files  were  dimensioned  according  to  longitude  (73),  latitude 
(11),  and  universal  time  (24).  Thus  Un  (  0,  0,  0)  corresponds  to  the  zonal  neutral  wind 
component  at  ±180  degree  longitude,  37.5  degree  latitude,  and  0  UT.  Each  data  file  thus 
contained  19,272  data  points. 

The  data  analysis  and  plotting  was  accomplished  almost  exclusively  by  the  IDL 
utility.  Figures  3.1  through  3.7  display  polar  plots  of  the  neutral  and  ion  winds  and  the 
forcing  terms  for  0,  6,  12,  and  18  UT.  The  relative  contributions  of  each  of  the  forces  is 
apparent  by  inspection  from  these  figures.  As  was  the  case  in  figure  2. 1 ,  the  ion  drag  and 
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pressure  gradient  forces  are  clearly  the  most  important,  while  the  viscous  drag  is  of  lesser 
significance. 

3.3  The  VSH  Model 

The  VSH  model  [Killeen,  et  al.,  1987]  was  developed  to  provide  a  more  “portable” 
version  of  the  TGCM.  The  VSH  is  essentially  a  decomposition  of  the  history  file  of  a 
specific  TGCM  run.  It  breaks  the  history  file  down  into  a  set  of  spherical  harmonic 
coefficients.  These  coefficients  can  then  be  applied  to  the  appropriate  polynomial 
expansion  to  return  the  value  of  the  desired  data.  There  are  two  features  of  the  VSH  which 
make  it  especially  attractive  to  this  study.  First  and  foremost  is  the  fact  that  it  is  extremely 
easy  to  use,  as  one  need  only  specify  a  short  set  of  parameters  in  a  FORTRAN  subroutine 
call  statement.  The  parameters  include  the  latitude,  longitude,  altitude,  time,  and  other 
geophysical  parameters.  This  flexibility  gives  the  VSH  its  second  attractive  feature,  as  it 
makes  it  very  easy  to  interpolate  between  the  grid  points  (both  spatial  and  temporal).  Thus, 
simulations  involving  a  satellite  orbiting  the  earth  in  virtually  any  orbit  are  very  easy  to 
perform.  The  output  of  the  VSH  includes  many  (but  not  yet  all)  of  the  outputs  of  the 
TGCM.  A  polar  plot  of  the  neutral  wind  as  computed  by  the  VSH  (including  a  short  list  of 
the  geophysical  parameters  used  as  input).is  shown  in  figure  3.8.  The  inputs  were  chosen 
to  duplicate  the  working  data  base  as  closely  as  possible.  A  comparison  of  this  figure  to 
figure  3.1  shows  that  the  VSH  output  is  acceptably  close  to  the  original  TGCM  data. 

3.5  Dynamics  Explorer  (DE-2)  Data  Base 

The  SPRL  library  at  the  University  of  Michigan  contains  an  on-line  data  base  of 
DE-2  measurements  for  over  8,000  orbits  of  the  satellite.  The  variety  of  data  available  is 
summarized  in  the  examples  contained  in  Appendix  D.  The  data  coverage  is  far  from  being 
continuous,  since  much  of  the  data  was  measured  in  situ  and  could  be  only  be  taken  during 
the  perigee  of  the  orbit.  In  addition,  not  all  the  satellite’s  instruments  were  activated  at  all 
times.  It  is  quite  often  the  case  that  a  particular  measurement  for  a  specific  location  at  a 
specific  time  is  not  available.  To  facilitate  orbit  selection  there  is  an  on-line  catalog 
available  which  lists  the  instrument  status  for  each  orbit.  For  this  particular  study,  the  only 
real  orbit  selection  criteria  were  that  the  satellite  provide  continuous  coverage  for  the  same 
latitude  band  on  three  consecutive  orbits.  Despite  this  apparently  lean  specification 
requirement,  the  data  base  could  provide  only  one  set  of  three  such  orbits. 

The  availability  of  the  DE  data  base  makes  possible  direct  comparisons  of  model 
output  to  satellite  observations.  Figure  3.9  displays  the  neutral  winds  along  two  segments 
of  DE  orbit  7217  along  with  simulated  TGCM  output.  The  satellite's  parameterized  orbit 


was  used  as  an  input  to  the  model  along  with  two  different  levels  of  geomagnetic  activity. 
As  is  evident  from  the  figure,  the  model  outputs  using  the  storm  level  inputs  show 
generally  good  agreement  with  the  satellite  data.  Thus  we  can  be  reasonably  assured  that 
the  TGCM  provides  a  realistic  picture  of  the  evolving  atmosphere. 
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Figure  3.9  DE  satellite  observed  neutral  winds  compared  with  TGCM  predicted  winds 
using  two  different  levels  of  geomagnetic  activity  as  input  Storm  input  case  agrees 
well  with  observed  winds. 
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Chapter  4  Experiments  with  the  Momentum  Equation 


The  chapter  contains  the  bulk  of  the  actual  work  in  this  project.  It  involves  using 
the  available  tools  to  try  to  experimentally  verify  the  momentum  equation  at  thermospheric 
altitudes.  The  first  section  we  consider  the  case  of  an  observer  fixed  at  some  specified 
location.  With  the  working  data  base  and  archived  data  available,  we  compare  the 
measured  acceleration  of  the  winds  and  compare  it  to  the  simultaneous  forces  acting  there. 
In  section  4.2  we  do  much  the  same  thing  for  the  case  of  a  special  class  of  moving 
observers  by  simulating  the  observations  of  polar  orbiting  satellites.  Finally,  in  section  4.3 
we  perform  the  same  analysis  on  a  set  of  archived  DE2  satellite  data. 


4.1  Fixed  Observer  Simulation 

In  this  section  we  seek  to  experimentally  verify  the  momentum  equation  using 
computer  generated  data  from  the  TIGCM.  The  momentum  equation  can  be  represented  in 
schematic  form  as  follows: 


f  =  SF„-Vr.  Vu 


Vv 


(4-1) 


Where  £  Fx  and  £  Fy  represent  the  sum  of  all  real  or  apparent  forces  acting  at  a  particular 
location.  These  forces  include  all  those  that  were  first  mentioned  in  the  second  chapter, 
namely,  ion  drag,  pressure  gradient,  Coriolis,  and  viscous  drag.  The  second  term  on  the 
right  hand  side  in  (4-1)  is  the  advection  term  where  Vr  is  the  relative  wind  and  is  defined  as 
the  difference  between  the  observer's  velocity  and  the  wind  velocity  (see  Apendix  A).  Our 
first  experiments  involve  fixed  locations,  which  is  equivalent  to  a  stationary  observer,  so 
the  value  of  Vr  is  simply  the  local  wind  velocity.  Quite  often  the  summation  of  the 
individual  force  terms  is  combined  with  the  advection  term  and  is  simply  referred  to  as  the 
total  force.  Appendix  A  explains  that  the  advection  term  is  not  a  real  force  in  the  same 
sense  that  the  other  terms  are,  but  is  necessary  to  account  for  relative  motion  between  the 
observer  and  the  medium  being  measured. 

If  all  the  forces  acting  on  the  neutral  atmosphere  are  accounted  for  in  (4-1),  then  we 
have  a  true  equality  between  the  observed  local  acceleration  (the  left  hand  side)  and  the  total 


force  term  (right  hand  side).  We  note,  however,  that  we  have  equality  only  in  the 
infinitesimal  limit  of  space  and  time.  Thus,  if  we  wish  to  experimentally  verify  the 
momentum  equation,  we  must  realize  that  the  best  we  can  hope  for  is  an  approximate 
equality.  We  will  proceed  in  a  fashion  similar  to  that  described  in  Appendix  A  and 
approximate  the  left  hand  side  of  (4.1)  by  calculating  the  slope  of  the  wind  component 
evolution  curves  for  a  fixed  location  and  comparing  this  value  with  the  total  force  computed 
by  the  model. 

Figure  4.1  depicts  a  set  of  microfiche-based  TGCM  data.  The  top  two  diagrams 
describe  a  24  hour  evolution  of  the  neutral  wind  components  at  the  location  shown  (42.5,  - 
85)  corresponding  to  the  approximate  location  of  Ann  Arbor.  This  data  set  was  obtained 
for  the  Z=1  pressure  level,  which  corresponds  to  approximately  300  km  in  altitude.  We 
observe  from  these  curves  that  the  neutral  wind  evolves  smoothly  over  time.  We  should 
therefore  expect  a  smooth  curve  drawn  through  the  discrete  data  points  to  be  a  good 
representation  of  a  continuous  set  of  data.  We  should  also  expect  to  be  able  to  calculate 
very  reasonable  approximations  to  the  time  derivative  of  the  neutral  wind  components. 

The  bottom  two  diagrams  show  the  time  evolution  of  all  the  forces  acting  at  this 
location.  The  breakdown  of  the  curves  is  shown  at  the  top  of  each  diagram  and  reads  as 
follows:  A  =  Ion  Drag;  B  =  Coriolis;  C  =  Pressure  Gradient;  D  *  Viscous  Drag;  and  E  = 
Advection.  There  is  also  a  curve  labelled  F  and  conveniently  represents  the  sum  of  the  A 
thru  E  curves.  We  observe  from  this  pair  of  diagrams  is  that  the  pressure  gradient  curve 
(C)  and  the  ion  drag  curve  (A)  are  clearly  the  dominating  forces.  We  also  note  that  these 
forces  tend  to  act  in  opposition  to  each  other,  so  that  the  total  force  curve  (F)  is  usually 
small  in  magnitude  and  does  not  venture  far  from  the  zero  force  line. 

The  values  for  the  neutral  wind  and  total  force  curves  can  be  taken  directly  from  the 
diagrams  and  the  results  are  shown  in  figure  4.2(a,c).  Data  points  were  taken  for  every 
hour  and  the  points  connected  using  a  natural  cubic  spline  function.  The  observed 
acceleration  of  the  wind  components  (figure  4.2(b))  was  taken  to  be  the  slope  of  the  wind 
evolution  curves  and  was  computed  using  interpolated  values  from  the  same  spline 
function.  The  momentum  equation  implies  that  the  total  force  curve  for  each  component 
should  be  very  similar  to  its  corresponding  observed  acceleration  curve. 

A  qualitative  comparison  of  figures  4.2(b)  and  4.2(c)  shows  that  the  curves  are 
indeed  very  similar.  We  can  get  a  feel  for  their  degree  of  similarity  by  integrating  the  total 
force  curve  with  respect  to  time  and  comparing  the  result  to  the  wind  evolution  curves.  To 
this  end  we  can  use  a  simple  Simpson's  rule  numerical  integration  on  an  odd  number  of  the 
discrete  values  for  the  total  force.  The  accuracy  of  this  method  can  be  tested  by  first 
performing  it  on  the  values  of  the  observed  acceleration.  The  results  in  figure  4.3(a)  show 
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that  the  integrated  values  (symbols)  fall  exactly  on  the  wind  component  curves  (smooth 
lines)  .  This  shouldn't  be  surprising,  since  Simpson's  rule  is  exact  for  cubic  functions, 
which  is  basically  what  a  spline  function  is. 

The  result  of  integrating  the  total  force  curve  is  shown  in  figure  4.3(b);  as  before, 
the  symbols  are  integrated  values  and  the  curves  are  the  actual  wind  c x  nponents.  While  in 
this  case  the  points  don't  lie  exactly  on  their  respective  curves,  the  general  trends  are 
acceptably  close  to  the  acceleration  curve.  The  fact  that  there  isn't  exact  agreement  could 
(in  this  particular  case)  be  the  result  of  errors  in  manually  digitizing  the  data  points,  errors 
in  using  the  cubic  spline  to  interpolate  the  actual  wind  and  forcing  evolutions,  or  even 
inaccuracies  in  the  graphs  of  the  archived  data.  The  numerical  integration  procedure 
shouldn't  be  considered  as  a  major  source  of  discrepancy,  since  it  was  shown  to  be  exact 
for  the  cubic  spline. 

This  experiment  can  be  repeated  for  virtually  any  location  using  the  working  data 
base  generated  by  the  TIGCM,  but  one  natural  location  to  consider  is  the  one  used  with  the 
archived  data.  A  summary  of  the  evolution  of  the  individual  forcing  terms  for  the  zonal  and 
meridional  neutral  wind  components  is  shown  in  figure  4.4.  The  breakdown  of  the  curves 
is  exactly  the  same  as  it  was  in  the  case  of  the  archived  data.  We  note  once  again  that  the 
pressure  gradient  and  ion  drag  terms  dominate  and  have  a  pronounced  tendency  to  oppose 
each  other.  Summing  of  all  the  individual  forces  produces  a  net  force  curve  (6)  that  is 
small  in  comparison  with  either  the  ion  drag  or  pressure  gradient  force. 

The  evolution  of  the  wind  components  for  the  Z=2  pressure  level  above  the  (42.5,  - 
85)  grid  point.is  shown  in  figure  4.5(a).  We  observe  once  again  the  very  smooth  nature  of 
the  changing  wind  components  and  claim  that  a  cubic  spline  fit  to  the  data  will  provide  a 
reasonably  accurate  approximation  to  a  continuous  data  set.  The  derivative  calculations  are 
again  made  using  interpolated  values  and  the  results  are  shown  in  figure  4.5(b),  and  they 
attest  to  the  relatively  inactive  nature  of  the  atmosphere  at  this  location.  Accelerations 
change  smoothly  and  the  magnitudes  are  small  (on  the  order  of  .005  m/sec2).  Based  on  the 
observed  accelerations  we  should  expect  the  total  force  curve  to  exhibit  similar 
characteristics;  namely,  small  magnitudes  varying  smoothly  with  time. 

The  total  forcing  curve  was  reproduced  in  figure  4.5(c)  using  the  same  scale  as  that 
of  figure  4.5(b).  It  is  immediately  apparent  that  the  two  figures  don’t  match  to  the  same 
degree  as  those  in  the  archived  case.  A  careful  comparison  of  the  two  figures;  however, 
shows  that  the  zonal  and  meridional  forcing  curves  actually  agree  quite  well  with  the 
observed  acceleration  curves  between  4  and  12  UT  (-2  and  6  SLT).  Beyond  that  time 
range  the  zonal  forcing  curve  grows  considerably  until  16  UT,  where  it  is  several  times  the 
magnitude  of  the  observed  acceleration.  It  then  drops  off  rapidly  until  it  shows  a  large 
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acceleration  in  the  opposite  direction.  Since  we  are  dealing  with  a  diumally  reproducible 
data  set,  the  net  effect  of  the  total  force  over  the  24  hour  period  ideally  should  be  zero.  The 
large  overshoot  of  the  total  force  after  12  UT  would  therefore  be  accompanied  by  a  large 
undershoot  prior  to  4  UT,  (where  the  total  force  curve  is  reasonable).  So  in  this  loose 
sense  the  total  force  curve  is  consistent,  but  it  doesn’t  explain  or  account  for  the 
discrepancies;  the  two  curves  should  be  similar. 

To  check  whether  this  discrepancy  is  a  general  condition  or  merely  a  local 
phenomenon,  the  analysis  was  repeated  for  three  latitudes  along  the  zero  degree  meridian. 
The  results  are  in  figures  4.6  through  4.8,  and  they  indicate  the  same  sort  of  response. 
During  the  eight  hour  period  prior  to  6  SLT  the  forcing  curves  for  both  components  are  in 
good  agreement  with  the  observed  acceleration;  the  meridional  component  generally  being 
the  better  of  the  two.  Beyond  6  SLT  the  forcing  curve  tends  to  show  greater  acceleration 
than  is  observed.  The  conclusion  is  that  the  apparent  imbalance  between  the  measured 
forces  and  the  observed  acceleration  is  not  limited  to  a  single  location. 

We  can  begin  to  get  a  feel  for  what  the  source  of  the  discrepancy  is  by  looking 
again  at  figure  4.4  and  noting  that  the  total  force  curve  (6)  rises  and  falls  in  concert  with  the 
ion  drag  curve  (1)  both  in  the  overshoot  of  region  excessive  force  (6  to  13  SLT)  and  in  the 
undershoot  region  (14  to  -2  SLT).  It  is  the  period  from  6  to  -2  SLT  when  the  ion  drag 
force  becomes  significant.  It  may  therefore  be  reasonable  to  assume  a  problem  with  the 
measurement  of  the  ion  drag  contribution  to  the  total  force.  To  explore  this  possibility,  it  is 
instructive  to  define  a  new  value  dubbed  the  "residual."  The  residual  is  the  value  the  ion 
drag  force  would  have  to  assume  to  provide  a  perfect  force  balance.  Expressed  numerically 
it  is  the  difference  between  the  observed  acceleration  and  the  sum  of  all  the  other  forcing 
terms  as  follows: 


Resx  =  -^-(lFx-IDx) 

Resy  =  lr_(lFy“IDy) 


(4-2) 


It  is  clear  from  these  equations  that  if  the  residual  were  used  in  place  of  the  ion  drag  force, 
then  the  total  force  and  local  acceleration  curves  will  exactly  coincide.  Plotting  the  residual 
and  ion  drag  forces  together  gives  immediate  indication  as  to  how  far  off  the  ion  drag  term 
really  is  (assuming  that  is  the  problem). 

The  ion  drag  and  the  residual  curves  are  plotted  together  in  figures  4.9  through 
4.12.  Where  the  residual  and  ion  drag  curves  are  close  together,  the  total  force  is  in  close 
balance  with  the  observed  acceleration.  When  the  two  curves  separate,  the  imbalance 
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grows.  A  quick  comparison  between  these  plots  and  those  of  the  ion  density  in  figure  4. 13 
shows  that  the  regions  of  force  imbalance  strongly  coincide  with  the  periods  of  increased 
ion  density,  which  is  one  of  the  determiners  of  the  strength  of  the  ion  drag  force. 

The  result  for  the  case  depicted  in  figure  4.9  (42.5,-85)  shows  that  the  zonal  ion 
drag  force  is  about  56%  greater  than  what  it  should  be  at  both  the  overshoot  and 
undershoot  extremes.  In  fact,  this  zonal  ion  drag  curve  appears  more  or  less  proportional 
to  the  residual  curve.  The  meridional  case,  on  the  other  hand,  shows  that  the  ion  drag  is 
actually  smaller  in  magnitude  than  the  residual.  This  means  that  if  the  ion  drag  force  were 
slightly  stronger  (about  20%),  then  the  meridional  total  force  would  balance  reasonably 
well  with  the  observed  acceleration. 

Examination  of  the  same  figures  for  the  zero  degree  meridian  (figs  4.10  -  4.12) 
reveals  a  similar  story.  Both  the  zonal  ion  drag  term  and  its  residual  exhibit  a  sinusoidal 
variation  between  6  and  0  UT.  The  zonal  ion  drag  appears  too  strong  at  the  extremes  of 
positive  and  negative  acceleration,  and  appears  roughly  proportional  to  the  residual  value 
especially  in  the  overshoot  region  where  ion  drag  is  greater.  It  would  also  appear  that  the 
phase  of  the  ion  drag  consistently  leads  the  residual  by  slightly  less  than  one  hour,  though 
it  is  unknown  what  significance  if  any  this  may  have. 

In  contrast  to  the  zonal  ion  drag's  sinusoidal  pattern  the  meridional  ion  drag  exhibits 
a  quasi-symmetric  behavior  about  local  noon.  The  figures  show  that  the  ion  drag  exceeds 
the  residual  in  magnitude  between  6  and  15  UT  with  maximum  discrepancy  occurring  near 
11  UT,  where  ion  drag  is  about  30%  too  strong.  That  the  residual  and  ion  drag  curves 
converge  at  15  UT  is  somewhat  of  a  surprise,  since  figure  4.13  shows  that  at  15  UT  the 
ion  density  is  still  more  than  half  its  daytime  maximum.  Another  interesting  feature  is 
found  in  the  curves  for  47.5  degree  latitude  (figure  4.12).  In  this  case  the  ion  drag  and 
residual  curves  appear  to  agree  slightly  better  than  for  the  lower  latitudes,  whereas  the  the 
zonal  drag  case  clearly  got  worse  with  increasing  latitude. 

To  try  explain  why  the  ion  drag  force  (as  calculated  by  the  diagnostics  processor)  is 
anything  other  than  what  it  should  be  (there  should  be  no  residual);  we  have  to  consider  a 
subtlety  in  the  TIGCM  model.  We  recall  that  the  archived  data  case  exhibited  good  force- 
acceleration  balance  for  the  entire  period  and  no  problem  with  ion  drag  was  evident. 
However,  the  particular  version  of  the  TGCM  that  created  that  data  used  previously 
archived  values  for  the  ion  drag  tensors  used  in  the  momentum  equation.  Since  the  same 
tensor  values  used  in  generating  the  global  wind  field  were  used  to  recover  the  ion  drag 
forces,  there  was  little  observed  discrepancy.  The  later  version  of  the  TIGCM  computed  its 
own  values  of  the  ion  drag  tensors.  The  tensor  values  recovered  by  the  diagnostics 
processor  could  conceivably  differ  slightly  from  those  used  during  the  model  run.  If  this  is 
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the  case,  then  the  ion  drag  force  thus  computed  will  be  different  enough  to  produce  a 
discrepancy  in  the  force  balance  equation  as  was  noted.  How  or  why  the  ion  drag  tensors 
should  differ  between  the  model  run  and  diagnostics  package  is  an  elusive  question  to  date. 

We  have  shown  in  this  section  that  the  momentum  equation  can  in  principle  be 
experimentally  verified  for  a  fixed  location  when  the  total  force  curve  is  known  to  be 
accurate.  We  have  also  shown  that  the  TIGCM  data  has  definite  but  limited  value  in 
demonstrating  this  balance.  We  must,  therefore,  exercise  caution  in  using  calculations 
from  this  data  where  the  ion  drag  force  is  dominant.  When  the  magnitude  of  the  ion  drag 
force  is  on  the  same  order  as  the  rest  of  the  terms  in  the  momentum  equation,  the  observed 
acceleration  tends  to  agree  with  the  total  force.  We  will  therefore  try  to  restrict  future 
computations  to  periods  when  the  ion  drag  term  is  at  a  minimum.  The  next  section 
describes  similar  experiments  for  the  case  of  a  moving  observer. 

4.2  Simulated  Satellite  Observations 

The  previous  section  considered  a  stationary  observer  at  a  fixed  location  in  space. 
The  results  showed  that  when  the  ion  drag  force  is  small,  the  agreement  between  the  local 
time  derivative  and  the  total  forces  is  actually  quite  good.  Problems  begin  when  ion  drag 
force  becomes  dominant  as  the  ion  density  grows  rapidly  after  local  sunrise.  This  section 
repeats  the  experiment  for  the  case  of  a  special  moving  observer:  a  satellite  orbiting  the 
earth. 

The  basic  problem  of  computing  the  total  derivative  of  a  measured  quantity  for  a 
moving  observer  is  discussed  in  great  detail  in  Appendix  A.  The  problem  states  that  for  an 
observer  moving  with  a  known  velocity  Vm,  he  will  observe  a  time  rate  of  change  in  that 
measured  quantity  according  to  the  following  equation: 

dQ  dQb 

-S^-dT-V  *Q  {Vr*Vw-Vm}  (4-3) 

The  first  term  on  the  left  hand  side  represents  the  rate  at  which  the  measured  quantity  is 
changing  within  the  medium  itself.  This  is  often  also  described  as  the  changes  observed 
while  following  the  motion  of  the  air.  It  is  instructive  to  think  of  this  term  as  the  rate  that  of 
change  within  the  airstream  at  the  observer's  location  and  at  a  specific  time.  To  avoid  any 
unnecessary  confusion  regarding  different  frames  of  reference  we  will  restrict  this 
discussion  to  a  single  reference  frame  (coordinate  system);  namely,  the  spherical  coordinate 
system  of  the  earth. 
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The  second  term  on  the  right  hand  side  (including  the  minus  sign)  is  the  advective 
term.  Its  function  is  to  relate  observed  changes  to  the  actual  changes  within  the  air  that 
result  from  physical  processes.  Note  that  the  advective  term  includes  both  the  wind 
velocity  (Vw)  and  the  observer's  velocity  (Vm),  and  both  are  measured  in  terms  of  the  well 


defined  coordinate  system  of  the  earth.  The  spatial  gradient  of  Q  is  likewise  computed 
within  the  earth's  coordinate  system. 

The  measurable  quantities  of  interest  in  this  study  are,  of  course,  the  neutral  wind 
components  u  and  .  The  physical  processes  responsible  for  changes  to  these  quantities 
within  the  medium  (the  airstream)  are  merely  the  forcing  terms  appearing  in  the  momentum 
equation.  The  form  of  the  momentum  equation  appropriate  for  the  earth's  rotating 
reference  frame  was  presented  in  chapter  2.  Replacing  the  first  term  on  the  right  hand  side 


of  (4-3)  with  these  forcing  terms,  we  come  up  with  the  following  equations 
du  1  3P  .  *1  32u  ,  r  .  .  uvtan  $ 
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(4-4b) 

Since  the  goal  of  this  section  is  to  demonstrate  force  balance  from  satellite  observations,  we 
must  consider  only  the  terms  that  are  calculable  from  satellite  observations.  The  first  term 
on  the  right  hand  side  of  equations  (4-4)  involves  the  pressure  gradient  force.  We  recall 
from  the  previous  section  that  the  pressure  gradient  force  was  largely  responsible  for 
balancing  the  ion  drag  force,  and  thus  it  is  an  important  contributor  to  the  total  force.  We 
also  note  that  to  compute  a  pressure  gradient  (or  any  gradient)  requires  simultaneous 
measurements  that  are  separated  spatially.  In  the  case  of  a  polar  orbiting  satellite  (such  as 
Dynamics  Explorer  2),  it  is  impossible  to  compute  a  zonal  pressure  gradient,  since 
measurements  are  only  available  along  the  track  of  the  satellite.  It  is  therefore  unlikely  that 
equation  (4-4a)  can  be  useful  in  demonstrating  force  balance,  so  calculations  involving  the 
acceleration  of  the  zonal  component  (u)  of  the  neutral  wind  will  not  be  considered. 

The  major  forcing  terms  of  the  meridional  component  ()  equation  (4-4b),  on  the 
other  hand,  can  be  calculated  (or  closely  approximated)  from  satellite  observations.  The 
pressure  gradient  term  in  (4-4b)  is  calculated  along  a  line  of  constant  longitude.  Typical 
polar  orbiting  satellites  travel  at  speeds  of  about  8  km/sec  along  a  line  of  nearly  constant 
longitude.  It  should  be  possible  in  principle  to  measure  the  pressure  along  the  satellite  track 
and  compute  the  meridional  pressure  gradient  from  these  measurements.  Of  course,  the 
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measurements  won't  actually  be  simultaneous,  but  the  actual  time  involved  is  small 
compared  to  local  accelerations  and  we  would  expect  this  to  be  a  reasonable  approximation 
to  the  actual  pressure  gradient. 

There  are  still  some  remaining  terms  which  cannot  be  calculated  from  satellite 
observations.  The  viscous  drag  term  contains  a  vertical  spatial  (second)  derivative  and 
cannot  be  computed.  The  viscous  drag  term  was  shown  in  chapter  3  to  be  a  very  small 
contributor  to  the  total  force.  Failing  to  include  this  term  should  therefore  have  minimal 
impact  on  momentum  equation  calculations.  The  advection  term  in  (4-4b)  can  be  expanded 
to  reveal  another  non-calculable  term: 

-(Vw"Vm)-  Vv  =  -[(u-  Vx)73^09X  +  (v_Vy)7|^]  (4-5) 


In  the  above  expression,  the  quantities  Vx  and  Vy  are  the  eastward  and  northward 
components  of  the  satellite's  velocity  with  respect  to  the  earth's  rotating  coordinate  system. 
For  a  satellite  whose  orbital  inclination  is  exactly  90  degrees,  Vy  is  merely  the  satellite's 
orbital  speed.  The  value  of  Vx  is  the  earth  corotational  speed  at  the  radius  of  the  orbit.  The 
first  term  in  the  brackets  on  the  right  hand  side  is  the  zonal  advection  of  the  meridional 
component  and  clearly  cannot  be  calculated.  The  second  term  is  the  meridional  advection 
of  the  meridional  term  (MAv)  and  can  be  approximated  in  much  the  same  manner  as  the 
pressure  gradient  force.  Whereas  the  viscous  drag  term  could  be  confidently  ignored,  it 
remains  to  be  seen  what  the  effect  of  neglecting  the  zonal  advection  term  will  be.  Yet  in 
spite  of  these  few  non-calculable  terms,  the  meridional  component  equation  would  seem  to 
be  well-suited  for  this  experiment. 

Collecting  all  the  terms  that  can  be  calculated  yields  the  working  form  of  the 
meridional  momentum  equation: 
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This  expression  would  at  first  appear  to  be  restricted  to  meridional  motion  only,  but  it  can 
be  suitably  adapted  to  represent  the  approximate  changes  observed  to  a  westward  moving 
observer  (satellite)  as  well. 

Consider  the  diagram  depicted  in  figure  4.13.  This  diagram  depicts  three 
successive  orbits  of  a  satellite  with  a  90  degree  inclined  orbit  and  a  90  minute  orbital 
period.  The  satellite  will  move  northward  along  lines  of  nearly  constant  longitude  during 
any  short  time  interval  of  a  given  orbit,  yet  progresses  22.5  degrees  (15  deg/hr  x  1.5  hrs) 
westward  after  one  complete  orbit.  Thus,  in  addition  to  the  obvious  case  of  a  northbound 
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observer  taking  continuous  observations,  there  is  an  analogous  case  of  a  westward  moving 
observer  recording  observations  every  90  minutes.  To  apply  this  analogy,  though,  it  will 
be  necessary  to  modify  (4-6)  by  removing  the  northward  component  of  the  satellite's 
velocity.  We  thus  now  have  the  following  two  forms  of  the  meridional  equation  to  work 


with: 
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The  first  expression  describes  the  observed  rate  of  change  of  component  v  along  a  line  of 
constant  latitude,  while  the  second  refers  to  the  rate  of  change  along  a  meridian.  It  cannot 
be  emphasized  enough  that  the  velocity  components  appearing  in  these  expressions  are 
measured  relative  to  the  earth's  coordinate  system  and  not  to  the  observer.  Thus,  in  the  last 
term  of  (4-7a)  v  is  the  wind's  actual  meridional  component.  So  with  two  working 
equations  in  hand,  it  is  now  possible  to  define  the  experiment. 


4.2.1  Simulation  1  Constant  Longitude 

This  experiment  examines  both  of  the  working  equations  by  simulating  an  orbiting 
satellite  collecting  data  along  its  track.  The  first  simulation  involves  a  satellite  moving 
along  a  meridian  recording  quasi-continuous  data.  The  second  simulation  tests  the  constant 
latitude  momentum  equation  by  looking  at  data  from  three  successive  orbits. 

The  first  experiment  demonstrates  the  effects  of  trying  to  apply  the  constant 
longitude  equation  to  a  short  orbit  segment.  This  is  accomplished  by  parameterizing  an 
ideal  orbit  along  a  single  meridian.  For  this  simulation  the  hypothetical  satellite  proceeds 
northbound  along  the  zero  degree  meridian  and  records  data  every  .5  degrees  of  latitude 
between  37.5  and  87.5  degrees.  This  corresponds  to  a  satellite  with  a  96  minute  orbital 
period  and  recording  data  every  8  seconds,  which  closely  approximates  the  performance  of 
the  DE-2  satellite.  Data  was  recorded  by  performing  time  interpolations  on  the  working 
TIGCM  data  base  at  5  degree  latitude  increments  from  37.5  to  87.5  degrees.  A  cubic 
spline  routine  was  then  used  to  interpolate  between  these  knot  points  at  .5  degree  intervals. 
Simulated  data  was  thus  obtained  for  Un,  Vn,  Vdrag,  Py,  VCor,  Geom  (curvature),  and 
MAv.  The  orbit  parameterization  and  recorded  data  are  summarized  in  table  C.  1  (see 
Appendix  C). 
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A  graph  of  the  recored  values  of  Vn  for  the  second  orbit  is  shown  in  figure  4.15(a). 
Each  symbol  represents  a  single  data  point  along  the  satellite's  track  as  data  is  recorded 
every  8  seconds.  The  graph  shows  that  it  took  the  satellite  only  800  seconds  to  travel  from 
37.5  to  87.5  degrees  in  latitude,  so  this  trace  should  be  a  close  approximation  to  the 
"snapshot"  of  a  meridional  distribution  of  Vn  along  the  zero  degree  meridian.  This  is 
illustrated  in  figure  4.15(b),  which  shows  the  distribution  of  Vn  along  the  zero  degree 
meridian  at  the  three  hourly  times  indicated.  Clearly,  the  profiles  don't  change  drastically 
from  one  hour  to  the  next,  so  no  significant  change  should  occur  along  the  meridian  during 
one  pass.  While  this  fact  makes  meridional  advection  calculations  possible,  it  has  dire 
consequences  as  far  as  the  momentum  equation  experiment  is  concerned. 

The  total  derivative  term  in  the  constant  longitude  equation  is  calculated  by 
computing  the  slope  along  the  Vn  trace  curve  at  each  point  of  interest.  This  is  easily  done 
using  finite  centered  differences,  and  the  results  are  shown  in  figure  4.16(a).  This  result  is 
then  compared  to  the  summary  of  the  momentum  equation  terms  shown  in  figure  4.16(b). 
The  most  striking  feature  of  these  figures  is  that  the  total  time  derivative  curve  from  4.16(a) 
is  virtually  identical  to  the  meridional  advection  term  of  4. 16(b).  The  reason  for  this  is  that 
the  satellite's  orbital  speed  of  some  8,000  m/sec  is  figured  in  the  momentum  advection 
calculation,  thus  resulting  in  a  dominating  term.  But  this  is  just  being  consistent  with  what 
was  stated  earlier;  namely,  that  the  meridional  wind  profile  isn't  going  to  change  very  much 
in  the  800  seconds  it  takes  for  the  satellite  to  complete  its  pass.  This  being  the  case,  it  is 
doubtful  that  the  forcing  effects  of  the  acceleration  on  the  airstream  can  be  detected  even  if 
the  known  effects  of  the  advection  are  subtracted  out.  In  practise,  there  will  be  a  certain 
amount  of  uncertainty  involved  with  the  wind  measurements,  which  will  carry  over  into  the 
advection  calculations.  It  is  likely  that  the  uncertainty  in  the  advection  term  alone  will  be 
greater  that  the  magnitude  of  the  instantaneous  acceleration  of  the  airstream.  More  on  the 
effects  of  uncertainty  is  discussed  in  the  next  section.  In  summary,  it  would  appear  that  the 
meridional  component  equation  isn't  well-suited  for  force  balance  analysis. 

4.2.2  Simulation  2  Constant  Latitude 

The  next  phase  of  the  experiment  involves  the  constant  latitude  momentum 
equation.  Recall  that  this  equation  is  analogous  to  a  westbound  observer  recording  wind 
data  every  90  minutes.  There  are  two  inherent  limitations  in  trying  to  experimentally  verify 
the  constant  latitude  momentum  equation.  As  was  previously  mentioned,  there  is  no 
practical  way  to  measure  the  zonal  gradient  of  Vn;  hence,  the  zonal  advection  term  cannot 
be  included.  In  addition,  observations  of  wind  velocity  are  necessarily  limited  to  one  per 
orbit  for  any  given  latitude.  The  total  time  derivative  of  the  wind  speed  must  be 
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approximated  using  a  centered  difference  with  a  time  step  on  the  order  of  100  minutes.  By 
making  appropriate  use  of  the  working  data  base  and  VSH  model,  it  is  possible  to  estimate 
the  likely  consequences  of  these  limitations. 

The  latter  limitation  is  usually  classified  as  a  limited  sampling  rate  problem,  and  it 
can  introduce  adverse  effects  in  at  least  two  ways.  Obviously,  the  sampling  rate  must  be 
sufficient  to  detect  changes  in  the  quantity  being  measured.  This  is  achieved  if  the 
sampling  rate  frequency  is  several  times  the  frequency  of  the  time  dependent  quantity  being 
measured.  In  addition  to  the  sampling  rate  limitation,  there  are  going  to  be  uncertainties  in 
measurements.  If  these  measurements  will  be  used  to  compute  derivatives  (as  in  this  case), 
then  the  uncertainties  will  carry  over  into  the  calculations  of  the  derivatives.  Both  these 
effects  can  be  demonstrated  rather  easily  with  the  VSH  model. 

To  demonstrate  the  sampling  rate  limitation,  a  simulated  observational  platform  was 
allowed  to  "orbit"  the  earth  in  a  westerly  direction  along  a  line  of  constant  latitude  and  with 
a  circumnavigational  period  of  24  hours.  The  starting  point  was  the  180  degree  meridian  at 
0  UT  at  an  altitude  of  395  km.  The  observer's  longitude  changed  at  a  rate  of  15  degrees 
per  hour  for  this  single  excursion  around  the  earth.  The  VSH  subroutine  was  called  at  130 
equally  spaced  points  along  the  track  with  the  appropriate  UT  for  each  point.  This  was 
done  for  several  different  latitudes  two  of  which  are  shown  in  figures  4.17(a)  and  4.18(a). 
The  continous  line  in  the  diagrams  represents  the  130  observations  of  Vn,  and  the 
diamonds  highlight  the  data  points  that  are  1.5  hours  apart.  By  inspection  of  these  figures 
it  would  appear  that  a  smooth  curve  drawn  through  the  highlighted  points  would  closely 
resemble  the  continuous  trace  curve.  The  sampling  rate  frequency  of  1.5  hours  is 
apparently  sufficient  to  detect  changes  in  the  winds.  The  next  step  is  to  determine  how 
accurate  an  estimate  to  the  total  time  derivative  can  be  achieved  with  these  data  points. 

The  error  bars  on  figures  4.17(a)  and  4.18(a)  describe  a  simple  5  percent  relative 
error  in  the  measurement  of  Vn.  This  value  was  arbitrarily  chosen  as  a  worst  case 
uncertainty  to  examine  the  effects  and  doesn't  represent  an  actual  estimate  of  the  error  likely 
to  be  encountered  in  the  model  values  or  in  satellite  measurements.  The  computations  of 
the  total  derivative  (dVn/dt)  are  shown  in  figures  4. 17(b)  and  4.18(b).  Again,  the  solid  line 
is  the  "actual"  time  derivative  taken  from  the  continuous  data  values,  and  the  diamonds  are 
the  approximations  to  the  derivative  computed  from  the  sampled  points.  The  error  bars  in 
this  case  result  from  the  5  percent  relative  error  assumption  for  the  neutral  wind 
measurement.  It  is  clear  from  these  figures  that  the  sampling  rate  limitation  precludes  exact 
measurement  of  the  total  time  derivative,  but  the  resulting  uncertainties  are  not  so  large  as  to 
preclude  a  reasonable  estimate  of  the  observed  total  time  derivative.  The  errors  from  the 


limited  sampling  rate  are  on  the  same  order  as  those  arising  from  the  uncertainties  in  the 
measurements. 

The  effect  of  neglecting  the  zonal  advection  and  viscous  drag  terms  can  be 
demonstrated  on  a  trial  basis  using  the  working  data  base.  This  was  done  by  selecting  a 
fixed  location  and  comparing  the  total  forces  present  during  a  24  hour  period  with  the  result 
obtained  by  subtracting  the  viscous  drag  force  and  zonal  advection.  Figures  4. 19  and  4.20 
show  these  results  for  two  locations.  It  should  be  noted  here  that  even  though  fixed 
locations  were  selected,  the  zoral  advection  appropriate  for  a  westward  drifting  satellite 
was  used.  In  this  comparison,  little  was  to  be  gained  by  allowing  the  location  to  change 
with  time. 

In  each  case,  the  solid  line  represents  the  total  forces  present  at  the  specified 
location  at  each  time.  The  dashed  lines  are  the  result  of  subtracting  the  zonal  advection  and 
viscous  drag  terms.  The  situation  presented  in  figure  4.19(a)  is  for  a  high  latitude  case  and 
shows  the  effect  of  large  variations  in  the  zonal  wind.  The  summary  of  the  zonal  wind 
component  for  this  location  is  shown  in  figure  4.19(b),  which  also  depicts  the  satellite  drift 
velocity  (Vmx).  Since  the  zonal  advection  term  of  (4-5)  involves  the  difference  in  the  zonal 
wind  and  drift  velocity,  the  effect  of  this  term  will  become  less  as  the  zonal  wind  speed 
approaches  that  of  the  satellite.  From  figure  4.19(b)  it  can  be  seen  that  this  happens 
between  10  and  13  UT.  It  is  during  this  same  period  that  the  difference  between  the  two 
curves  from  figure  4.19(a)  is  small. 

The  first  four  hours  in  this  case  represent  the  other  extreme,  where  the  zonal  wind 
adds  to  the  effect  of  advection.  This  enhancement  effect  is  reflected  in  figure  4.19(b). 
where  there  is  a  noticeable  maximum  in  the  difference  between  the  two  curves.  This  same 
effect  is  demonstrated  in  figures  4.20  for  the  case  of  a  lower  latitude  and  weaker  winds.  In 
this  case  the  effect  of  the  zonal  wind,  though  evident,  is  not  as  significant  because  the  zonal 
wind  magnitude  is  much  smaller  than  the  drift  velocity.  It  is  also  evident  from  the  last  six 
hours  of  the  low  latitude  case  that  the  wind  speed  is  not  the  only  determining  factor.  Here 
the  difference  between  the  two  curves  is  at  a  minimum  for  the  24  hour  period,  yet  the  zonal 
wind  is  not  significantly  different  than  in  earlier  periods.  The  magnitude  of  the  gradient  of 
the  meridional  component  must  be  contributing  as  well,  but  as  was  previously  mennoned, 
cannot  be  measured  from  satellite  data.  The  conclusion  to  draw  from  these  figures  is  that 
the  zonal  advection  term  can  be  effectively  neglected  when  the  zonal  wind  speed  is  close  to 
the  drift  velocity  of  the  satellite.  When  working  with  actual  satellite  data,  it  is  not  always 
possible  to  be  very  selective  about  which  locations  to  consider. 

The  procedure  for  this  phase  of  the  experiment  was  to  parametrize  the  orbits 
depicted  figure  4.14  and  use  appropriate  interpolations  between  the  grid  points  of  the 
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working  data  base  to  create  a  set  of  observations  for  this  hypothetical  satellite.  The 
parametrization  was  arbitrarily  chosen  so  that  the  satellite  passed  through  (42. 5N,  OW)  at 
precisely  6  UT,  and  the  90  minute  orbital  period  defined  the  rest  of  the  coordinate  pairs. 
The  reference  starting  point  on  each  orbital  pass  was  taken  to  be  30  degrees  north,  and  data 
was  recorded  every  8  seconds  after  the  satellite's  latitude  first  exceeded  37.5  degrees  (the 
southern  limit  of  the  working  data  base).  Data  values  were  linearly  interpolated  from  the 
data  base  at  latitudes  37.5,  42.5,  47.5,  52.5,  57.5,  and  62.5;  the  values  for  longitude  and 
UT  being  consistent  with  the  orbit  parameterization.  These  data  values  were  then  inserted 
into  a  cubic  spline  routine  to  interpolate  data  values  for  each  of  the  latitudes  where  data  was 
recorded.  The  orbital  parameterization  along  with  retrieved  values  of  Un,  Vn,  Vdrag, 
Vcor,  Py,  Geom  (curvature),  and  MAv  are  shown  in  table  C.2. 

The  measurements  of  Vn  from  these  three  orbits  are  shown  in  figure  4.21(a).  The 
dotted  line  in  this  figure  is  the  middle  orbit  and  is  the  one  along  which  the  momentum 
calculations  are  made.  Note  that  the  poleward  gradient  of  Vn  increases  with  time.  This  is 
to  be  expected,  since  the  starting  time  is  close  to  6  UT,  and  we  recall  from  the  fixed 
observer  case  that  the  ion  drag  force  increased  substantially  after  6  UT.  It  is  tdso  evident 
from  the  figure  that  the  time  r^e  of  change  of  Vn  increases  substantially  with  latitude  as 
time  goes  on.  This  effect  is  not  seen  at  low  latitudes,  where  the  latest  orbit  actually  showed 
weaker  meridional  winds  during  the  third  orbit. 

The  sum  total  of  the  measurable  forces  (and  advection)  are  shown  in  figure  4.21(b). 
Also  shown  in  this  figure  is  the  curve  representing  the  total  time  derivative  of  Vn.  Ideally, 
the  solid  curve  and  dashed  curves  should  be  close  together  because  each  curve  represents 
opposite  sides  to  the  constant  latitude  momentum  equation.  The  solid  line  is  the  left  hand 
side  and  is  computed  at  each  latitude  by  taking  centered  differences  in  the  meridional  wind. 
The  dashed  curve  is  the  right  hand  side  and  is  the  sum  of  all  the  available  force  terms 
including  advection.  Since  the  two  curves  bear  only  superficial  resemblance,  it  might  be 
proposed  th.  t  the  zonal  advection  of  the  meridional  wind  is  significant  here.  To  explore 
that  possibility  a  third  curve  was  added  to  the  figure,  and  is  the  sum  of  the  available  forces 
and  the  zonal  advection  term.  Again,  the  two  curves  bear  only  superficial  resemblance  and 
defy  explanation. 

The  values  of  the  measurable  forces  during  the  second  orbit  are  plotted  together  in 
figure  4.22(a).  The  importance  of  the  pressure  gradient  and  ion  drag  force  terms  is  clearly 
evident.  From  the  earlier  work  with  the  fixed  location  momentum  equation,  it  was  noted 
that  the  ion  drag  and  pressure  gradient  forces  tend  to  oppose  each  other.  From  the  figure  it 
would  at  first  appear  that  this  is  not  quite  the  case,  as  both  forces  tend  to  grow  in  the 


positive  sense  north  of  50  degrees  latitude.  However,  if  the  entire  zero  degree  meridian  is 
shown  (as  in  figure  4.22(b)),  the  opposing  nature  of  the  forces  is  evident. 

As  in  the  case  of  the  fixed  observer,  verifying  the  momentum  equation 
experimentally  remains  an  elusive  goal.  At  least  in  the  former  case  it  was  shown  to  work  in 
the  absence  of  strong  ion  drag  force.  In  this  case  it  doesn't  seem  to  verify  at  all.  In  the 
present  case,  the  total  forces  present  at  40  degrees  latitude,  for  example,  would  predict  an 
instantaneous  acceleration  of  the  wind  of  about  .002  m/sec^.  In  actual  fact,  the  wind 
showed  a  slight  deceleration  between  the  first  and  third  orbits.  Even  if  advection  is  fully 
accounted  for,  the  total  forces  tend  to  exceed  the  predicted  acceleration.  An  unusual  twist 
to  the  present  case  is  that  it  is  the  pressure  gradient  force  and  not  the  ion  drag  that  is  causing 
the  apparent  problem,  since  the  total  force  is  skewed  in  the  positive  sense  relative  to  the 
observed  acceleration.  The  last  experiment  involves  actual  satellite  data;  so  if  there  is  a 
problem  with  using  model- generated  data  to  perform  momentum  equation  experiments, 
then  perhaps  the  satellite  data  could  yield  more  favorable  results. 

4.3  Momentum  Equation  Calculations  from  Satellite  Data 

In  the  preceding  sections  the  momentum  equation  calculations  were  based  entirely 
on  model-generated  data.  In  these  section  we  perform  essentially  the  same  computations 
with  the  data  taken  from  the  DE-2  data  base.  The  first  problem  one  encounters  when  trying 
to  do  any  experiment  with  this  data  is  to  find  a  satisfactory  set  of  orbits.  By  satisfactory  we 
mean  orbits  that  contain  enough  of  the  data  you  need,  in  the  right  locations  in  the  orbit  for 
your  application,  and  contains  enough  suitable  orbits  to  make  time  studies  possible.  For 
this  experiment  we  required  zonal  neutral  winds,  meridional  neutral  and  ion  winds,  ion  and 
neutral  temperatures,  electron,  atomic  oxygen  and  molecular  nitrogen  densities.  To  avoid 
problems  with  high  latitude  forcing,  it  was  preferable  to  have  data  at  low  to  mid  latitudes. 
Finally,  to  make  time  derivative  computations  as  accurately  as  possible,  we  required  three 
consecutive  orbits'  worth  of  this  data.  Almost  miraculously,  this  stringent  set  of 
requirements  was  filled. 

Out  of  a  data  base  consisting  of  thousands  of  orbits,  there  was  only  one  set  of  three 
that  contained  full  data  coverage  over  a  desirable  latitude  band.  Orbits  8025  through  8027 
filled  the  requirements  perfectly.  They  contained  all  the  required  data,  where  it  was 
needed,  and  for  as  long  as  it  was  needed.  Figure  4.23  shows  a  schematic  diagram  of  the 
track  of  these  three  orbits  over  the  south  geographic  pole.  The  vectors  along  the  satellite 
track  represent  the  neutral  wind.  Note  that  is  representation  of  the  south  pole  is  as  viewed 
from  the  north  pole.  Thus,  east  is  in  a  counterclockwise  direction  (the  direction  of 
increasing  local  time).  Note  also  that  the  neutral  wind  vectors  show  a  strong  westerly 


34 


component  in  the  55  to  70  degree  latitude  range.  As  was  previously  mentioned,  this  would 
tend  to  reduce  the  effect  of  zonal  advection  and  improve  the  results. 

The  procedure  for  this  experiment  is  very  similar  to  that  of  the  previous  section. 
The  lower  portion  of  figure  4.23  shows  the  track  of  the  three  orbits  between  latitudes  -40 
and  -60.  As  before,  the  change  in  the  meridional  wind  from  orbit  8025  to  8027  is 
compared  to  the  total  forces  acting  along  the  track  of  orbit  8026.  The  only  real  difference 
between  this  experiment  and  the  simulation  (aside  from  no  longer  needing  to  parameterize 
the  orbits)  is  that  all  the  forces  must  be  computed  from  the  satellite  data  rather  than 
interpolated  from  the  working  data  base.  A  summary  table  of  all  the  DE-2  satellite  data 
within  this  20  degree  latitude  band  of  interest  as  well  as  plots  of  the  data  listed  in  the  first 
paragraph  can  be  found  in  Appendix  D.  Descriptions  of  the  force  computations  and  plots 
of  the  forces  along  the  track  will  be  found  in  the  paragraphs  to  follow. 

4.3.1  Ion  Drag  Force 

The  ion  drag  force  is  calculated  using  the  following  formula  taken  from  the 
momentum  equation: 

Vdrag  =  vn.(v.-v„)  =  (4-8) 

Equation  (4-8)  looks  deceptively  simple,  but  this  apparent  simplicity  belies  the  complex 
nature  of  computing  the  momentum  transfer  collision  frequency,  vin.  The  ion-neutral 
collision  frequency  is  actually  a  combination  of  the  collision  frequencies  between  the 
oxygen  ion  and  the  two  dominant  neutral  species  listed  above.  The  formulas  for  these  two 
collision  frequencies  are: 

^iN2  =  ^3  n^2  (4-9a) 


V 


iO 


=  C  ]  nO  T1/2)  1  -  C2Iog  ]0Tr)2 


(4- 9b) 


where: 


ViN2 

ion  -  nitrogen  collision  frequency 

vio 

ion  -  oxygen  collision  frequency 

nN2  = 

nitrogen  number  density 

nO 

oxygen  number  density 

ne 

election  (ion)  density 

Tr  = 

(Ti  +  Tn)/2 

Cl  =  Resonant  collision  frequency  coefficient 

C2  =  Oxygen  collision  frequency  constant 

C3  =  Non-resonant  collision  frequency  coefficient 

In  the  above  equations  the  constants  Cl,  C2,  and  C3  are  taken  from  Banks  and  Kockarts 
[1973]  and  are  1.7  x  3.67  x  10'^,  0.064,  and  6.82  x  IQ'10,  respectively.  A  detailed 
analysis  on  the  effect  of  uncertainties  in  the  measured  quantities  is  presented  in  Appendix 
B.  The  result  of  the  error  analysis  shows  that  if  a  1  percent  relative  error  is  assumed  for 
each  of  the  quantities  going  into  the  computation  of  the  ion  drag,  the  resulting  error  in  the 
calculation  is  less  than  10  percent.  This  would  probably  represent  a  best  case  estimate  of 
the  ion  drag. 


4.3.2  Pressure  Gradient  Force 

The  pressure  gradient  force  is  in  principle  very  easy  to  calculate.  At  each  point 
along  the  track  of  the  satellite  measurements  are  taken  of  the  neutral  densities  and  the 
neutral  temperature.  Thus  all  the  ingredients  are  there  to  calculate  a  pressure  along  each 
1/2  degree  of  latitude  as  long  as  data  is  being  recorded.  If  the  satellite  moved  along  at  a 
constant  altitude,  then  the  pressure  gradient  force  could  be  calculated  with  the  following 
formula: 


P 


y 


(4-10) 


where  p  is  the  mass  density  of  the  gas.  In  practise  the  satellite’s  altitude  usually  changes  as 
it  moves  along  its  track.  To  correct  for  the  changing  altitude  we  use  the  barostatic  relation: 

P cor  =  P(z)exp  [  -  (z  -  z  0)/H  ]  (4-11) 


where:  Par 

P(z) 
A) 

H 


Pressure  corrected  to  reference  altitude 
Pressure  at  satellite  altitude 
Reference  altitude 
Mean  scale  height  along  track 


The  mean  altitude  along  the  track  between  -40  and  -60  was  taken  to  be  the  reference 
altitude.  After  correcting  the  pressure  values,  a  single  application  of  a  3-point  smoothing 
routine  was  applied  to  dampen  small  scale  variations  between  successive  data  points.  The 
smoothed  corrected  pressure  along  the  track  of  orbit  8026  is  shown  in  figure  4.24(a).  The 
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figure  shows  a  definite  trend  toward  higher  pressure  as  one  moves  poleward,  but  there  are 
also  shorter  scale  undulations  in  pressure  along  the  track. 

Computing  the  pressure  gradient  involves  little  more  than  simply  performing  finite 
difference  methods  along  the  track,  but  there  is  a  potential  problem  involved  when 
measurement  uncertainties  are  involved.  If  we  assume  as  before  that  measurements  are 
accurate  to  a  relative  error  of  1  percent,  then  the  error  involved  in  computing  the  pressure  is 
shown  by  the  error  bars  in  figure  4.24(a).  The  error  appears  to  be  small,  but  close 
examination  reveals  the  error  to  exceed  the  amount  of  pressure  change  between  adjacent 
data  points.  This  condition  has  severe  consequences  when  it  comes  to  using  centered 
differencing  in  computing  the  pressure  gradient,  as  can  be  seen  in  figures  4.24(b).  The 
error  bars  in  this  figure  result  from  the  uncertainties  in  the  pressure  values  and  the  distance 
between  the  data  points.  The  latter  quantity  appears  in  the  denominator  of  the  expression 
for  the  absolute  error  of  the  pressure  gradient  formula  (see  Appendix  B).  One  could  halve 
the  calculation  error  if  instead  of  using  every  point  to  compute  the  pressure  gradient,  every 
other  point  were  used  instead.  There  is  a  practical  limit,  of  course,  to  how  many  points 
you  can  skip  and  still  have  a  representative  derivative  calculation. 

The  actual  uncertainty  in  this  particular  pressure  gradient  calculation  is  probably 
considerably  less  than  what  is  shown  in  figure  4.24(b).  One  of  the  advantages  in  having  a 
nearly  continuous  set  of  data  points  is  that  it  allows  you  to  perform  smoothing  operations 
and  thus  improve  upon  the  uncertainty  of  individual  data  points.  Thus,  the  errors  shown  in 
figure  4.24  may  tend  to  be  on  the  high  side. 

4.3.3  Time  Derivative  Computation 

Computing  the  total  time  derivative  of  the  meridional  wind  component  is  basically 
as  simple  as  differencing  the  wind  measurements  from  orbits  8027  and  8025  and  dividing 
by  the  time  difference.  That,  along  with  some  enhancements,  is  just  what  was  done.  To 
make  sure  that  wind  values  from  identical  latitudes  were  differenced,  the  raw  wind  data 
was  input  into  a  cubic  spline  routine  and  interpolated  values  for  each  1/2  degree  latitude 
increment  were  obtained.  Three  point  smoothing  was  performed  on  the  resulting  wind 
values  and  the  differencing  calculations  were  performed.  The  values  of  the  interpolated 
winds  are  shown  in  table  4-1,  as  are  the  force  calculations  along  the  second  orbit.  A  plot  of 
the  time  derivative  is  shown  in  figure  4.25(a).  The  error  bars  once  again  result  from  the 
assumption  of  a  1  percent  relative  error  in  the  velocity  measurements. 
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4.3.4  Force  Comparison 

A  comparison  plot  of  the  calculated  forces  is  shown  in  figure  4.25(b).  The 
uncertainty  in  the  pressure  gradient  force  calculation  dominates  the  picture  as  a  whole.  The 
ion  drag  force  curve  appears  to  be  well-behaved,  which  lends  itself  to  a  certain  amount  of 
credibility.  The  undulations  in  the  pressure  gradient  force  are  not  reflected  in  the  ion  drag 
force,  which  is  what  one  would  come  to  expect  based  on  previous  results.  There  is, 
however,  an  encouraging  rise  in  the  pressure  gradient  force  at  about  -57  degrees  which 
roughly  corresponds  to  an  opposite  move  in  the  ion  drag  force.  It  would  have  been  better, 
of  course,  had  it  been  more  closely  in  phase  to  its  ion  drag  counterpart,  especially  in  light 
of  the  comparison  of  the  total  force  to  the  observed  acceleration  of  the  wind. 

4.3.5  Force  Balance  Analysis 

The  somewhat  disappointing  results  of  this  experiment  are  summarized  in  figure 
4.25(c).  The  scale  of  the  graph,  so  chosen  to  keep  the  total  force  curve  within  its 
boundaries,  effectively  kills  the  curve  for  the  observed  acceleration.  Actually,  the 
acceleration  curve  wrongly  suggests  that  nothing  was  happening  with  the  neutral  wind 
during  this  period  when,  in  fact,  it  had  gone  up  and  then  down  between  orbits  8025  and 
8026,  nearly  returning  to  the  same  values  and  yielding  the  small  values  for  acceleration. 
The  obvious  conclusion  from  this  experiment  is  that  no  conclusion  is  possible,  at  least  not 
within  the  level  of  detail  presented  here. 
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Figure  4.1  Archived  TGCM  data  showing  diurnal  evolution  of  the  zonal(a)  and 
meridional  (b)  neutral  wind,  and  the  zonal  (c)  and  meridional  (d)  force 
components. 
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Figure  4.2(a)  Evolution  of  neutral  wind 
components  as  appearing  in  archived 
TGCMdata. 


Figure  4.2(b)  Observed  acceleration  of 
neutral  winds  as  calculated  from 
curves  above. 


Figure  4.2(c)  Total  force  dau  taken  directly 
from  archived  TGCM  data.  The 
curves  should  resemble  those  of  fig 
4.2(b)  above. 
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Figure  4.4  Forcing  summary  for  same  location  as  in  archived  case,  but  using  data 
from  the  TIGCM  working  data  base. 
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Figure  4.5(a)  Diurnal  evouldon  of  neutral 
wind  as  taken  from  TIGCM  working 
data  base  ( 42.5,  -85 ). 
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Figure  4.5(b)  Measured  acceleration  of  the 
neutral  wind  depicted  above. 
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Figure  4.5(c)  Total  forces  acting  at  same 
location  as  above  figures. 


Figure  4.5 
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Figure  4.6(a)  Oiumal  evoultion  of  neutral 
wind  as  taken  from  T1GCM  working 
data  base  (  37.5, 0). 


Figure  4.6(b)  Measured  acceleration  of  the 
neutral  wind  depicted  above. 


Figure  4.6(c)  Total  forces  acting  at  same 
location  as  above  figures. 
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Figure  4.7(a)  Diurnal  evouldon  of  neutral 
wind  as  taken  from  TIGCM  working 
data  base  ( 42.5, 0). 
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Figure  4.7(b)  Measured  acceleration  of  the 
neutral  wind  depicted  above. 
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Figure  4.7(c)  Total  forces  acting  at  same 
location  as  above  figures. 


Figure  4.7 
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Figure  4.9  Measured  ion  drag  force  compared  with  the  ion  drag  that  would  provide  perfect 
balance  of  forces  (ie  the  “residual”)  for  location  ( -42.5,  -85  ). 
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Figure  4.8 
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Comparing  Meridional  Ion  Drag  and  Residual 


Figure  4.10  Measured  ion  drag  term  compared  with  the  residual  for  location  (  37.5,  0). 
Meridional  residual  is  greatest  during  the  sunlit  hours. 


48 


Acceleration  (m/sec2)  Acceleration  (m/sec2) 


Comnarincr  Zonal  Ion  Drac  with  Residual 


0  1  2  3  4  5  6  7  8  9  10  1 1  12  13  14  15  16  17  18  19  20  21  22  23 

Universal  Time 


Cnmnarinor  MariHinnal  Tr»n  Drau  anrl  Rasirlual 


0  l  2  3  4  5  6  7  8  9  10  1 1  12  13  14  15  16  17  18  19  20  21  22  23 

Universal  Time 


Figure  4.1 1  Ion  drag  and  residual  term  compared  at  ( 42.5, 0  )  shows  the  zonal  term  with 
increasing  residual  magnitude  after  local  noon. 
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Comparing  Zonal  Ion  Drag  with  Residual 
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Comparing  Meridional  Ion  Drag  and  Residual 


Figure  4.12  Ion  drag  and  residual  at  (  0,47.5  )  shows  the  ion  drag  has  at  least  twice  the 
magnitude  of  the  residual  term  in  the  late  afternoon.  Meridional  term  as  much 
closer  by  comparison. 
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Electron  (Ion)  Deneity 


4.13  Ion  density  temporal  evolution  for  locations  of  previous  four  figures.  Ion 
density  increase  corresponds  separation  of  residuals  from  ion  drag  force  in 
meridional  forcing  curves. 
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Figure  4.15  Measurements  of  Vn  along  the  zero  degree  meridian  from  a 
northbound  satellite  (a)  and  for  three  different  hourly  sets  (b). 
Measurements  along  track  approximate  the  instantaneous  profiles. 
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Figure  4.16  Calculation  of  the  instantaneous  rate  of  change  of  Vn  (a)  and  a 
composite  of  calculable  forces  along  satellite  track  (b).  The  satellite’s  orbital 
speed  causes  advection  term  to  dominate  and  closely  resemble  that  in  (a). 
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Figure  4.17  Vn  and  dVn/dt  for  a  hypothetical  westbound  observer.  Plotting 
symbols  are  drawn  1.5  hrs.  Solid  lines  in  (b)  are  exact  derivatives  while 
plotting  symbols  are  centered  difference  approximations.  Five  per  cent 
relative  error  bars  in  (a)  give  rise  to  computational  errors  in  (b). 
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Figure  4. 1 8  Same  depiction  of  Vn  and  dVn/dt  for  a  high  latitude.  Larger  variations 
in  Vn  lead  to  greater  accelerations.  Uncertainty  in  calculating  acceleration 
by  centered  differencing  is  on  the  same  order  as  the  uncertainty  resulting 
from  the  measurements. 
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Figure  4.19  Comparison  of  the  total  meridional  forcing  to  that  excluding  the 
viscous  drag  and  zonal  advection  (a).  Agreement  is  best  where  zonal  wind 
speed  approaches  the  satellite  westward  drift  velocity  (b). 
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Figure  4.20  Same  comparison  as  in  Figure  4. 19  showing  that  zonal  wind  speed  has 
a  noticeable  albeit  reduced  effect  when  zonal  winds  are  light 
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Figure  4.21  Results  of  constant  latitude  experiment  showing  Vn  along  the  three 
orbital  bands  (a).  Acceleration  between  the  first  and  third  orbits  shows  the 
effect  of  increasing  forcing  with  latitude  (b),  but  is  far  below  total  forcing 
curve.  Effect  of  restoring  zonal  advection  term  also  shown. 
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Figure  4.22  Summary  of  forcing  along  second  simulated  orbit  (a).  Pressure 
gradient  force  and  ion  drag  force  along  entire  meridian  (b)  show  familiar 
opposing  characteristics. 
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Figure  4.24  Smoothed  pressure  values  along  second  orbit  (a).  Plotting  symbols 
denote  observation  points.  Meridional  pressure  gradient  (b)  showing  error 
bars  resulting  from  assumed  uncertainties  in  pressure  calculation. 
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Figure  4.25  Summary  of  results  showing  comparison  of  observed  acceleration  (a) 
to  the  indivduaJ  forces  (b)  and  directly  to  the  total  (c).  Observed 
acce'eration  is  far  less  than  instantaneous  forces. 


Chapter  5  Conclusions  and  Remarks 


The  objective  of  this  thesis  was  to  demostrate  the  validity  of  the  momentum 
equation  based  on  both  theoretical  and  experimental  data.  To  this  end  I  obtained  a  large 
amount  of  theoretical  data  from  a  single  run  of  the  NCAR  general  circulation  model.  I  also 
had  at  my  disposal  volumes  of  microfiched  archive  data  from  several  previous  runs  of  an 
early  version  of  the  TGCM  as  well  as  vast  amounts  of  DE-2  satellite  data. 

My  approach  was  very  simple  in  principle.  The  momentum  equation  in  its  eularian 
form  states  that  the  instantaneous  acceleration  of  the  wind  at  some  location  is  exactly  equal 
to  the  sum  of  all  the  forces  that  are  acting  upon  the  airstream  at  that  same  location  for  that 
same  instant  in  time.  In  Appendix  A  I  carefully  derived  the  physical  relationship  between 
moving  and  stationary  observers  and  how  their  measurements  would  differ  in  a  very 
precise  way.  I  also  showed  here  that  it  was  possible  to  relate  the  slope  in  the  trace  of  a 
moving  observer's  chart  recorder  to  the  changing  values  of  the  measured  quantity  within 
the  medium  and  account  for  the  observer's  relative  motion.  This  set  the  stage  for  the 
experiment  to  include  measurements  and  computations  based  on  satellite  data. 

In  section  4.1  I  applied  the  predictions  of  the  momentum  equation  to  a  fixed 
location  in  space  and  used  microfiche  data  from  an  early  TGCM  run.  In  spite  of  the 
relatively  crude  manner  in  which  data  was  transferred  from  the  microfiche  the  data  analysis 
software,  the  results  for  this  phase  were  extremely  encouraging.  The  local  acceleration  did 
indeed  appear  to  be  equal  to  the  sum  of  the  terms  in  the  momentum  equation.  This 
enthusiasm  was  short-lived  as  I  proceeded  to  the  next  phase  of  the  experiment. 

When  I  tried  to  apply  the  same  analysis  to  the  working  data  base,  the  results  were 
not  nearly  as  good.  I  found  that  the  forces  began  to  exceed  the  measured  acceleration  as 
soon  as  the  ion  drag  force  became  appreciable.  When  the  ion  drag  force  was  small  in 
comparison  to  the  sum  of  the  other  forces,  the  momentum  equation  did  appear  to  v'erfy. 
The  apparent  discrepancy  could  not  be  readily  attributed  to  impulsive  ion  forcing  because 
the  drag  force  was  a  smoothly  varying  function  of  time  whose  functional  form  could  be 
approximated.  It  was  decided  to  proceed  with  the  next  phase  of  the  experiment  and  choose 
locations  where  ion  drag  was  expected  to  be  small. 

Section  4.2  involved  satellite  simulations  in  an  analogous  case  to  that  of  the 
stationary  observer.  This  time,  hov  ver,  I  intentionally  didn’t  include  the  forces  that  could 
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not  be  measured  from  a  real  satellite.  As  expected,  the  results  were  not  as  favorable  as 
those  of  the  fixed  observer  cases.  What  was  not  expected  was  the  degree  towhich  the 
momentum  equation  failed  to  verify.  Even  when  the  missing  terms  were  restored  to  the 
force  calculation,  the  results  were  marginal  at  best. 

Finally,  in  section  4.3  I  applied  the  analysis  to  real  satellite  data.  The  data  appeared 
to  be  very  reasonable  and  free  from  noise.  It  was  unusually  complete  and  smooth  curves 
of  raw  data  were  obtained.  However,  when  all  the  calculations  were  made,  the  forcing 
curves  beared  no  resemblence  whatever  to  the  acceleration  curve.  When  plotted  on  the 
same  scale,  the  total  forcing  curve  showed  significant  time-dependent  structure,  yet  no 
variations  in  the  accerlation  curve  could  be  discerned. 

In  spite  of  the  apparent  lack  of  success  in  this  experiment  I  demonstrated  that  most 
of  the  significant  forces  in  the  momentum  equation  could  be  calculated.  The  problem  when 
it  comes  to  trying  to  use  the  total  force  to  compute  local  accelerations  is  that  the  total  force 
us  usually  small  in  comparison  with  either  the  ion  drag  and  pressure  gradient  or  both.  The 
difference  between  the  two  is  often  very  small  and  it  is  the  difference  that  provides  the 
driving  force  for  the  atmosphere.  Any  uncertainties  in  the  computation  of  either  of  the 
forces  could  very  well  be  sufficient  to  make  local  acceleration  predictions  virtually 
meaningless. 

I  still  believe  it  is  possible  to  use  satellite  data  in  a  moving  observer  analogy; 
perhaps  the  momentum  equation  is  not  the  one  to  use  to  experiment  with  this.  The  energy 
equation  is  also  involved  in  the  TGCM  and  may  be  a  candidate  for  a  similar  study. 
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Appendix  A  The  Physical  Meaning  of  the  Total  Derivative 


In  making  applications  of  the  momentum  equation  one  has  to  be  absolutely  certain 
of  the  physical  meaning  of  all  the  terms  involved.  While  it  may  seem  a  rather  simple 
concept,  the  total  time  derivative  of  the  wind  components  has  historically  been  a  source  of 
considerable  confusion  or  misconception.  It  is  the  purpose  of  this  section  to  remove  any 
confusion  regarding  the  total  derivative  (as  it  is  used  throughout  this  thesis)  and  in  the 
process  attach  a  physical  meaning  to  this  term. 

It  is  safe  to  say  that  establishing  a  clearly  defined  coordinate  system  is  fundamental 
to  any  discussion  regarding  the  measurement  of  some  physical  quantity  at  different  times 
and  at  a  variety  of  locations.  For  the  present  discussion  consider  the  rectangular  coordinate 
system  depicted  in  figure  A.l.  The  axes  are  labelled  with  arbitrary  values,  but  it  may  be 
instructive  to  consider  the  units  to  be  statute  miles  and  the  grid  pattern  to  represent  an  area 
on  the  earth's  surface.  The  X  and  Y  directions  naturally  correspond  to  the  east  and  north 
directions  respectively.  With  the  coordinate  system  thus  defined,  we  can  calculate  all 
positions  and  velocities  with  respect  to  its  origin. 

Within  this  coordinate  system  we  will  be  recording  measurements  of  some  scalar 
quantity  Q  that  is  a  characteristic  of  the  air.  Q  can  be  any  continuous  measurable  quantity, 
but  for  our  purposes  it  is  useful  to  think  of  it  as  the  temperature  of  the  air.  The  labelled 
dashed  lines  in  figure  la  represent  an  analysis  of  the  spatial  distribution  of  Q  at  some  point 
in  time  we  have  designated  as  t=0.  From  figure  A.l  it  is  apparent  that  the  value  of  Q  is  a 
function  of  X  and  Y  with  a  local  maximum  at  location  (10,15).  Figures  A.  1(a)  and  A.  1(b) 
depict  the  same  coordinate  system  and  Q  distribution  at  the  two  future  times  indicated.  It  is 
clear  from  the  figures  that  the  Q  distribution  is  moving  eastward  as  time  progresses;  the 
physical  interpretation  being  that  the  airmass  is  moving  (with  respect  to  our  coordinate 
system)  and  carries  its  values  of  Q  with  it.  By  comparing  the  analyses  in  each  figure  it  also 
becomes  apparent  that  the  distribution  itself  is  changing  as  it  moves  eastward.  The  physical 
interpretation  of  this  is  that  the  airmass  itself  is  changing  in  the  quantity  Q  as  time  goes  on. 
If  Q  is  the  air  temperature,  then  we  can  think  of  the  airmass  as  a  whole  as  gradually 
warming  up  with  time.  Clearly  then,  the  value  of  Q  depends  not  only  on  the  location  of  the 
observer  within  the  coordinate  system  but  on  the  time  of  the  observation  as  well. 


Within  the  coordinate  system  are  three  observers;  one  observer  I  -d  at 

location  (30,20),  while  the  other  two  are  moving  within  the  coordinate  s>  sterr.  \  mobile 
observer  (M)  is  located  at  (90,30)  at  t=0  and  is  travelling  in  a  westerly  direction  I  he  third 
observer  (B)  is  riding  in  a  free  flying  balloon  and  as  such  is  constrained  to  follow,  the  wind; 
he  is  thus  stationary  in  relation  to  the  airmass.  Figure  A. 2  is  a  repeat  of  figure  A.  1  (b)  with 
the  paths  of  observers  M  and  B  superimposed  on  the  coordinate  system.  Note  that  at  time 
t=3  the  locations  of  all  three  observers  coincide. 

It  is  further  presumed  that  all  three  observers  are  taking  continuous  measurements 
of  Q  and  are  recording  them  on  a  chart- type  recorder  as  shown  in  figure  A. 3.  Note  that  at 
time  t=3  the  values  of  Q  recorded  by  each  observer  is  the  same  (same  location  and  time). 
Note  further  that  while  the  instantaneous  value  of  Q  is  the  same  for  all  three  observers  at 
this  point ,  the  slopes  of  the  three  curves  are  quite  different  at  this  point.  In  other  words, 
the  rate  of  change  of  Q  is  different  for  all  three  observers  even  though  the  value  itself  is  the 
same.  The  next  step  is  to  try  to  calculate  the  value  of  the  rate  of  change  for  all  three 
observers. 

Let  Q(x,y,t)  represent  the  spatial  and  temporal  dependency  of  the  quantity  Q.  Next 
consider  a  differential  change  in  Q  and  expand  it  according  to  the  rules  of  differential 
calculus  as  follows: 


dQ=^<lx  +  |Sdy+^dt 


at 


(A-l) 


Dividing  both  sides  of  this  expression  by  dt  yields: 

dQ  dQ  dx  dQ  dy  dQ 
dt  dx  dt  dy  dt  dt 


(A-2) 


The  left  hand  side  of  this  expression  is  just  the  time  rate  of  change  of  quantity  Q  as 
determined  by  an  observer.  Physically,  this  is  the  slope  of  the  trace  on  the  observer's  chart 
recorder.  We  recognize  the  quantities  dx/dt  and  dy/dt  on  the  right  hand  side  as  the  X  and  Y 
components  of  observer's  velocity  (as  defined  in  our  coordinate  system  ).  We  can  consider 
the  case  of  the  mobile  observer  (M)  and  rewrite  (A-2)  as  shown. 

dQm  dQ  dQ  dQ 

dt  ~  "air v™  +‘5y"Vmy+"ar  (a-3) 


The  subscript  denotes  values  that  are  unique  to  the  particular  observer.  The  last  term  on  the 
right  hand  side  is  the  time  rate  of  change  of  Q  the  observer  would  record  if  his  velocity 


68 


were  zero;  in  other  words,  the  slope  on  his  trace  if  he  were  fixed  at  some  location.  In  our 
example,  this  term  is  the  slope  of  observer  L's  chart  trace.  We  can  rewrite  (A-3)  in  vector 
form  by  recognizing  the  first  two  terms  on  the  right  hand  side  as  the  scalar  product  of  a 
velocity  and  spatial  gradient  of  Q. 


dQm 

dt 


=  V 


m 


VQ  + 


aQ 

■3T 


(A-4) 


This  is  the  familiar  definition  of  the  total  time  derivative  of  Q.  As  we  have  already  noted, 
its  physical  interpretation  can  be  thought  of  as  the  slope  of  the  trace  of  continuous  Q 
observations  taken  by  an  observer  who  is  moving  with  velocity  V  within  a  spatially 
dependent  field  of  Q.  This  value  is  a  combination  of  the  changes  in  Q  occurring  at  that 
particular  location  and  the  changes  he  observers  as  a  result  of  his  motion.  Clearly,  the  total 
derivative  is  unique  to  the  observer,  but  the  local  time  derivative  (the  second  term  in  the 
right  hand  side  of  (A-4) )  is  observer  independent.  We  can  come  up  with  a  quantitative 
expression  for  the  local  derivative  if  we  repeat  (A-3)  from  the  viewpoint  of  the  balloon- 
borne  observer  (B). 


dQt 

dt 


iQv  . 

dx 


3 

ay  %  +  at 


(A-5) 

Since  B  is  forced  to  follow  the  airmass,  we  can  substitute  the  components  of  the  wind 
velocity  into  (A-5)  and  come  up  with  the  following  equivalent  expression: 

dQ, 


dt 


aQ  aQ  aQ 

dx  u  +  ay v  +  at 


(A-6) 


Solving  equation  (A-6)  for  aQ/dt  yields  the  following  expression. 


dQ 

~3T 


dQ, 


dt 


aQ  aQ 

-£TU  +  57v 


(A-7) 


Or  equivalently  in  vector  form: 

dO  dQ. 

ir  =  -dr-vr'VQ  (vr-vw-vm}  (a-8) 

The  first  term  on  the  right  hand  side  (in  addition  to  being  the  slope  of  B's  chart  trace  in 
figure  A.3)  represents  the  rate  at  which  the  quantity  Q  is  changing  within  the  airmass  itself. 
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If  Q  is  the  temperature,  then  this  term  reflects  the  rate  at  which  the  airmass  (or  more 
correctly,  this  part  of  the  airmass)  is  warming  due  to  physical  heating  processes. 

The  vector  quantity  on  the  right  (including  the  minus  sign)  is  commonly  called  the 
advection  of  Q.  It  is  the  quantity  responsible  for  observed  changes  in,  say,  temperature 
that  result  purely  from  the  motion  of  the  airmass  past  this  fixed  location.  We  can  substitute 
expression  (A-7)  for  the  local  derivative  into  (A-3)  to  obtain  after  rearranging: 


dQm  _  dQb 

dt  dt 


dQ  dQ  1 

(u  -Vmx)-5T  +  (v  ~Vmy)-^- 


(A-9) 


Once  again  changing  to  the  equivalent  vector  form: 


3Q  <** b 

c >r  dt 


Vr-  VQ 


(A- 10) 


The  vector  Vr  is  defined  as  the  difference  between  the  observer’s  velocity  and  that  of  the 
wind  and  is  commonly  referred  to  as  the  relative  wind.  Note  if  the  observer's  velocity  is 
zero,  then  his  total  derivative  assumes  the  role  of  a  local  derivative  and  we  have  the 
situation  in  equation  (A-8).  If  on  the  other  hand  the  wind  velocity  is  zero,  then  B's  total 
derivative  becomes  a  local  derivative  and  we  are  back  to  equations  (A-3)  and  (A-4). 

Equation  (A- 10)  gives  us  a  form  of  the  total  derivative  that  is  very  useful  in 
interpreting  the  momentum  equation.  If  instead  of  Q  we  use,  say,  the  zonal  wind 
component,  then  the  first  term  on  the  right  hand  side  in  (A- 10)  becomes  dut/dt  and  is  the 
total  derivative  of  u  following  the  motion  of  the  air.  In  other  words,  it  is  acceleration  of  the 
air  at  the  point  of  the  observation  that  is  the  result  of  the  sum  of  all  real  forces  acting  at  that 
point.  These  forces  are  the  individual  terms  in  the  momentum  equation  that  do  not  include 
advection.  The  advection  term  is  the  second  term  on  the  right  hand  side  and  is  a  function  of 
both  the  local  wind  velocity  and  the  observer's  velocity  within  the  coordinate  system. 
Equation  (A- 10)  stated  in  words  unequivocally  assigns  the  physical  meaning  to  the  total 
derivative  as  follows: 

The  instantaneous  rate  of  change  of  a  component  of  the  wind  for  an  observer 
moving  within  a  well  defined  reference  frame  is  equal  to  the  total  acceleration 
of  that  component  of  the  airflow  at  that  location  minus  the  scalar  (dot) 
product  of  the  wind  velocity  relative  to  the  observer  and  the  locally  observed 
gradient  in  the  component  being  measured. 
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Figure  A.  I  Hypothetical  distribution  of  a  scalar  value  Q  at  three  different  time*.  Diagram 
show  location*  of  stationary  and  two  moving  observers  (L,  M,  and  B). 
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Figure  A.2  Situation  at  time  t  =  3  showing  all  observers  at  same  location. 
Trajectories  of  moving  observers  have  been  superimposed  on  the  diagram. 
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Figure  A.3  Superimposed  chart  traces  from  the  three  observers.  The  curves  cross  at 
t*3  because  all  observers  were  at  the  same  location.  Different  slopes  at  this 
point  result  from  differing  observer  velocities. 
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Appendix  B  Calculation  of  Errors  in  the  Momentum  Equation 


Throughout  this  project  little  emphasis  was  placed  on  estimating  the  actual  absolute 
error  in  the  raw  satellite  data.  Conversely,  much  was  done  to  examine  the  effects  of  errors 
in  the  subsequent  calculations  involving  these  measured  quantities.  Specifically,  we 
consider  the  errors  in  the  forcing  term  calculations  arising  from  an  estimated  relative  error 
in  the  raw  data.  The  treatment  of  error  propagation  in  calculations  is  treated  in  Greenberg 
[1975]  and  is  applied  in  this  section. 

In  a  function  involving  several  variables  the  error  in  the  calculation  arising  from 
uncertainties  in  the  variables  can  be  estimated  as  follows: 


F(X1>X2’X3"‘) 


8F  =  -4p-8x .  + 

dxt  1 


3F  x  .  9F  g 
5x.5x2+  3x,6x3 


(B-l) 


in  the  above  expression,  8F  represents  the  absolute  error  in  the  calculation  of  F  resulting 
from  the  absolute  uncertainties  in  its  constituent  quantities  (5xi,  8x2,  8x3,  ...).  This 
expression  actually  estimates  the  maximum  possible  error,  since  a  straight  forward  sum  of 
the  individual  uncertainties  is  used  to  make  the  estimate.  It  is  unlikely,  however,  that  in 
any  one  calculation  all  of  the  individual  quantities  would  simultaneously  be  in  error  by  the 
maximum  range  in  their  uncertainties.  For  this  reason  we  choose  to  use  the  probable  error 
in  these  calculations.  For  an  absolute  error  consisting  of  a  sum  of  several  individual  errors, 
then  the  probable  error  is  computed  as  follows: 

<5FW  a/£<K>2  (B-2) 


where  5xn  are  the  absolute  errors  in  the  measurement  of  quantity  x,> 

If  a  calculation  involves  th  product  of  several  variables,  each  containing  its  own 
uncertainty,  then  the  relative  error  in  the  value  of  the  function  is  the  sum  of  the  individual 
relative  errors.  The  probable  error  involving  such  a  product  is  computed  as  shown. 


(B-3) 
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The  applications  of  these  fundamental  principles  to  the  present  problem  are  addressed  in  the 
paragraphs  that  follow. 

All  momentum  calculations  in  this  project  are  based  on  data  that  is  available  from 
the  DE  satellite  data  base.  Among  this  includes  the  ion  temperature  (Tj ),  neutral 
temperature  (Tn),  atomic  oxygen  density  (nO),  nitrogen  density  (nN2),  electron  density 
(ne),  neutral  meridional  wind  (Vn),  ion  meridional  wind  (Vi),  and  neutral  zonal  wind  (Un). 
No  attempt  is  made  here  to  estimate  the  absolute  errors  likely  to  be  encountered  in  these 
data;  the  purpose  of  this  section  is  to  observe  how  momentum  calculations  are  affected  by 
errors.  Thus,  each  of  the  quantities  listed  is  assumed  to  have  an  inherent  relative  error  of 
.01  (1%). 


B.l  Observed  Acceleration  (AVn) 

A  fmite-centered  differencing  technique  is  used  to  estimate  the  value  of  the  time  deri  ative 
(or  acceleration)  of  the  observed  neutral  wind.  For  each  latitude  an  observation  is  made 
once  each  orbit.  Thus  if  we  record  wind  observations  (Vo,  Vj,  V2, ...)  every  At  seconds, 
the  acceleration  is  given  by: 

“•!T4(vrvJ  <B-4> 


Applying  (B-l)  to  this  function  yields  the  following  expressions  for  the  absolute  possible 
and  probable  errors: 


8AV„  = 


8AV„ 

- -SV 

av„  ov2i 


3AV 


-5V, 


3V0  ”'0| 


<«Av„)po!S  =  2ir(8v2)  +  1ir(Sv0) 

<*AVn)^=  ^(SV^-MSV,,)2] 


(B-5a) 

(B-5b) 

(B-5c) 


where  8V0  and  8V2  are  the  absolute  errors  taken  from  the  value  of  V2  and  Vo  and  the 
assumed  relative  error  of  1  percent. 


B.2  Coriolis  Force  (VCor) 

The  Coriolis  force  is  easily  calculated  from  the  measured  wind  as  follows: 

VCor  =  -  fU  n 


7  4 


(B-6) 


Where  f  is  the  Coriolis  parameter  and  depends  only  on  latitude.  From  this  formula  and  the 
above  discussion  the  absolute  error  in  the  Coriolis  term  becomes: 


SVCor  =  fSU  n 


where  5U  n 


V  un  J 


(B-7) 


B.3  Meridional  Advection  (MAv) 

The  meridional  term  of  the  momentum  advection  can  be  approximated  from  the  available 
DE  data: 

dVn 

MAvS-Vv-^T=-VnGradVn  (B-8) 


The  relative  error  in  the  above  calculation  of  advection  is  then: 


f  8MA  v  ^ 

2 

1 

(  SGradV  n  ^ 

T 

l  mav  J 

prob 

Li  Vn  J 

-r 

l  GradVn  J 

(B-9) 


The  first  term  in  brackets  is  the  relative  error  in  Vn  and  is  either  estimated  or  determined 
experimentally.  Vn  is  measured  along  the  satellite  track  and  is  virtually  continuous,  so 
finite-centered  differencing  can  be  used  to  compute  GradVn  as  follows: 

GradV"=  W(V2-Vo)  (B-10) 


The  probable  error  in  this  calculation  is  then: 


SGradV  n 


1 

2  Ay 


+ 


(B-ll) 


Because  the  Ay  term  appears  in  the  denominator,  one  must  use  care  in  computing  gradients 
along  the  track  of  the  satellite.  If  the  smallest  Ay  is  used,  the  resulting  uncertainty  in  me 
calculation  might  be  larger  than  the  expected  magnitude  of  the  gradient.  If  Ay  is  too  large, 
then  the  approximation  to  the  spatial  derivative  becomes  less  accurate. 


B.4  Pressure  Gradient  Force  (Py) 

The  pressure  gradient  force  is  calculated  with  the  following  formula: 
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P>  =  ‘f?(fr)  =-<P)_'<GradP) 


(B-12) 


The  subsequent  expression  for  the  error  is  found  by  applying  (B-l).  The  result  is: 


<SPy)poS, 


(apA 


p+ls&spJ5GradP 


(B-13a) 


(SP  )  =  Ol^P5p+  ‘scradP 

v  poss  pz  K 


(B-l  3b) 


(8Py)  =  GradP-^y  + 

prob  |_V  P 


2  2 

Spl  ( SGradP^ 

o2  l  p  J 


(B-l 3c) 


The  neutral  mass  density  (p)  and  uncertainty  are  readily  calculated  from  the  raw  data: 

p  =  (6.01  x  1020)  (rrO  •  nO  +  mN^nN^  (kg/m3)  (B-14a) 


-If  2 

8p  =  V6.01  x  1020)  [(mO  •  8nO)  +  (  mN  2  •  8nN  2)  "  (B-l 4b) 

where  mO  and  mN2  are  atomic  and  molecular  weights  of  oxygen  and  nitrogen.  The 
pressure  gradient  is  calculated  in  the  same  manner  as  the  meridional  velocity  gradient 
above.  The  resulting  erro  *  in  this  calculation  is: 


6GradP=  2^(5P2)2+,5Po)2] 


(B-l  5) 


As  before,  the  value  of  Ay  should  be  carefully  chosen  to  ensure  the  magnitude  of  the 
uncertainty  does  not  exceed  that  of  the  gradient.  The  pressure  and  its  uncertainty  are 
calculated  from  the  ideal  gas  law: 

P  —  n  n  k  ^Tn  (B- 1 6a) 


SP  _  Sn^  8T n 
P  “  nn  *  T 

11  n 


(B-16b) 


The  second  term  on  the  right  hand  side  of  (B-16b)  is  the  estimated  relative  error  of  the 
neutral  temperature.  The  first  term  on  the  right  hand  side  is  the  relative  error  in  the  neutral 
density: 

nn=nO+nN^  =>  8nn  =  j^(8nO)  +  (8nN2)  (B-17) 

Finally,  the  absolute  probable  error  in  the  pressure  calculation  is  given  by: 


(B- 18) 


B.5  Ion  Drag  (ID) 

The  ion  drag  force  and  resulting  possible  relative  error  can  be  calculated  with  the  following 


expressions: 


ID  =  v  V.  nr 

in  in  r 


8v .  SV .  8n 


ID  v 


Ji+  _iL+  __I 


Vin  " r  (B- 19) 


where:  Vjn  =  ion-neutral  momentum  transfer  collision  frequency 

Vin  =  (Vi-Vn) 

nr  =  ne/nn 

The  last  two  terms  on  the  right  hand  side  are  easily  calculated  as  follows: 


(B-20a) 


5v,„=[(«v,)2+(8Vr)2] 


(B-20b) 


The  error  in  the  collision  frequency  requires  many  more  computations,  the  total  ion-neutral 
collision  frequency  being  the  sum  of  the  collision  ion-oxygen  and  ion-nitrogen  collision 
frequencies: 


v  =  v,,  +  v,,  - 

m  O  iN  2 


r  2  2i 

8tln  =L(5v<>)^+  (  8V  iN  i) 


(B-21) 


The  formulas  for  calculating  the  ion-nitrogen  collision  frequency  and  uncertainty  are  fairly 
simple: 


(B-22) 


v  =  C  n 
IN2  L3  N2 


^ViN2“  C3^nN2 


The  resonant  ion-oxygen  collision  frequency  requires  a  more  complex  formula: 


vd  =  c1"oti1'2b2 


(B-23) 


where:  Tr  =  .5  (Tj  +  Tn) 

B  =  1  -  C2  Logio(Tr) 

Equation  (B-l)  must  be  rigorously  applied  to  this  expression  as  follows: 


5v°  ■  a^8"0  +  W. :ST' 


(B-24a) 


where  8Tr  =  .5[(8T,)2  +  (ST„)2] 


(B-24b) 


Performing  the  indicated  partial  differentiations  eventually  leads  to  the  following  expression 
for  the  absolute  error  in  the  ion-oxygen  collision  frequency: 


f  2  2s  1/2 

8vo=-  [cj/'VSno]  +  -^(B2-  1.740^81,  ■ 

.  ^  r 


(B-25) 


B.6  Total  Force  (TOTy) 

The  absolute  probable  error  for  the  total  force  is  found  by  summing  the  individual  errors  in 
the  usual  way: 

r  2  2  2  2  i1/2 

8T0Ty  =[  (8ID)  +  (8Py)  +  (5MAV)  +  (8VCor)  J  (B-26) 
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Appendix  C  Simulated  Satellite  Data 


This  section  contains  listings  of  the  simulated  data  and  calculated  quantities  in  the 
satellite  simulations  described  in  sections  4.2.1  and  4.2.2.  Because  of  their  size  the  data 
tables  were  included  in  an  appendix  rather  than  in  the  main  body.  Table  C.l  contains  the 
data  from  the  constant  longitude  portion  of  the  experiment.  Table  C.2  contains  the  data  and 
computations  from  the  constant  latitude  case. 
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SIMULATED  DATA  RECORDED  ST  A  RORTBBOUHC  SATELLITE  MOVI1G  ALOSG  EERO  DEGREE  NE1IDIA1 


■•(•) 


0 


at 

Ua 

Vn 

*7 

VCor  Vdias 

G«oa 

HAt 

6Ta/4t 

7.5 

-5.21 

35.9 

4. 624-03 

4.62E-04  -7. 644-03 

-3 

044-06 

-l 

734-02 

-1.734-0 

1.0 

-4.05 

35.7 

4.61E-03 

5.46E-04  -7.66E-03 

-4 

254-06 

-1 

734-02 

-1.741-0 

(.5 

-4.44 

35.6 

4.61E-03 

6.27E-04  -7.631-03 

-5 

554-06 

-1 

764-02 

-1.771-0 

9.0 

-7.41 

35.4 

4.61E-03 

7. 054-04  -7.591-03 

964—06 

-1 

411-02 

-1.011-0 

9.5 

-4.33 

35.3 

4.60E-03 

7. 794-04  -7.561-03 

494-06 

-1 

474-02 

-1.004-0 

0.0 

-9.01 

35.1 

4.S9E-03 

4. 504-04  -7.511-03 

014-05 

-1 

964-02 

-1.971-0 

0.5 

-9.45 

35.0 

4.S4E-03 

9.144-04  -7.461-03 

-1 

141-05 

•  2 

064-02 

-2.074-0 

1.0 

-10.2 

34.4 

4.S7E-03 

9.444-04  -7.401-03 

-1 

351-05 

-2 

191-02 

-2.201-0 

1.5 

-10.4 

34.6 

4. 551-03 

1.044-03  -7.341-03 

541-05 

-2 

351-02 

-2.361-0 

2.0 

-11.3 

34.4 

4.S3C-03 

1.104-03  -7 . 275-03 

-1 

721-05 

-2 

531-02 

-2.541-0 

2.5 

-11.4 

34.2 

4 .  S1E-03 

1.164-03  -7.201-03 

-1 

924-05 

-2 

721-02 

-2.741-0 

3.0 

-12.3 

34.0 

4.441-03 

1.224-03  -7.121-03 

-2 

114-05 

-2 

901-02 

-2.911-0 

3.5 

-12.4 

33.7 

4.4SE-03 

1.244-03  -7.041-03 

-2 

304-05 

•  3 

021-02 

-3.031-0 

4.0 

-13.2 

33.5 

4.42E-03 

1.334-03  -6.951-03 

-2 

494-05 

-3 

061—02 

-3.041-0 

4.5 

-13.5 

33.2 

4.39E-03 

1.341-03  -6.451-03 

•  2 

644-05 

-3 

051-02 

-3.061-0 

5.0 

-13.4 

33.0 

4.36E-03 

1.421-03  -6.751-03 

-2 

451-05 

-2 

961-02 

-2.971-0 

5.5 

-14.1 

32.4 

4.32E-03 

1.464-03  -6.631-03 

-3 

011-05 

-2 

401-02 

-2.421-0 

4.0 

-14.3 

32.5 

4.29E-03 

1.494-03  -6.511-03 

—  3 

141-05 

-2 

541-02 

-2.591-0 

6.5 

-14.4 

32.4 

4.261-03 

1.524-03  -6.341-03 

-3 

264-05 

-2 

241-02 

-2.291-0 

7.0 

-14.5 

32.2 

4.23E-03 

1.544-03  -6.241-03 

-3 

344-05 

-1 

914-02 

-1.921-0 

7.5 

-14.4 

32.0 

4.21E-03 

1.554-03  -6.091-03 

-3 

391-05 

-1 

464-02 

-1.471-0 

0.0 

-14.3 

31.9 

4 .194-03 

1.554-03  -5.931-03 

-3 

391-05 

-3 

664-03 

-9.701-0 

4.5 

-14.1 

31.9 

4.14E-03 

1.541-03  -5.761-03 

-3 

351-05 

-4 

251-03 

-4.271-0 

9.0 

-13.4 

31.9 

4.17S-03 

1.511-03  -5.591-03 

*3 

251-05 

531-03 

1.541-0 

9.5 

-13.3 

31.9 

1. 174-03 

1.444-03  -5.411-03 

•  3 

111-05 

771-03 

7.401-0 

0.0 

-12.4 

32.0 

4 . 141-03 

1.431-03  -5.241-03 

-2 

904-05 

451-02 

1.451-0 

0.5 

-12.1 

32.1 

4 . 19E-03 

1.361-03  -5 . 041-03 

-2 

641-05 

141-02 

2.191-0 

1.0 

-11.2 

32.4 

4.21E-03 

1.271-03  -4.931-03 

-2 

324-05 

901-02 

2.991-0 

1.5 

-10.2 

32.6 

4.24E-03 

1.171-03  -4.791-03 

-X 

964-05 

451-02 

3.461-0 

2.0 

-9.06 

33.0 

4.24E-03 

1.041-03  -4.671-03 

-1 

564-05 

401-02 

4.421-0 

2.5 

-7.47 

33.4 

4.32E-43 

4.451-04  -4.541-03 

-1 

131—05 

451-02 

5.441-0 

3.0 

-6.06 

33.9 

4.34E-03 

7.021-04  -4.521-03 

101-02 

7.131-0 

3.5 

-4.25 

34.5 

4. 454-03 

4.941-04  -4.491-03 

-3 

631-06 

621-02 

4.661-0 

4.0 

-2.26 

35.3 

4.S2E-03 

2.621-04  -4.441-03 

-i 

051—06 

10 

0.10 

4.5 

-0.112 

36.2 

4.60E-03 

7.441-06  -4.491-03 

-2 

631-49 

12 

0.12 

5.0 

2.19 

37.3 

4.69E-03 

-2.641-04  -4.531-03 

-1 

021-06 

15 

0.15 

5.5 

4.64 

34.7 

4.74C-03 

-5.631-04  -4.591-03 

—4 

651-06 

17 

0.14 

4.0 

7.22 

40.2 

4. 471-03 

-4.771-04  -4.661-03 

141-05 

21 

0.21 

4.5 

9.92 

42.1 

4.961-03 

-1.201-03  -4.761-03 

-2 

201—05 

24 

0.25 

7.0 

12.7 

44.2 

9. 051-03 

-1.551-03  -4.471-03 

-3 

691-05 

29 

0.29 

7.5 

15.6 

46.4 

9.141-03 

-1.921-03  -5.001-03 

-3 

691-05 

33 

0.33 

4.0 

14.6 

49.7 

9.221-03 

-2.301-03  -5.141-03 

-4 

241-05 

30 

0.30 

4.5 

21  .  7 

52.9 

9.291-03 

-2.694-03  -5.241-03 

-1 

141-04 

42 

0.42 

9.0 

24.4 

56.4 

9 . 331-03 

-3.091-03  -5.411-03 

-1 

511-04 

45 

0.45 

9.5 

27.4 

60.2 

9.341-93 

-3.501-03  -5.501-03 

-1 

951-04 

47 

0.44 

0.0 

30.9 

64.2 

9.301-43 

-3.901-03  -5.541-03 

-2 

451-04 

49 

0.50 

0.5 

33.9 

64 . 2 

9.20B-03 

-4.291-03  -5.311-03 

-3 

011-04 

51 

4.51 

1.0 

34 . 7 

72.4 

9.03E-03 

-4.671-03  -5.401-03 

-3 

611-04 

52 

0.52 

1.5 

39.5 

76.7 

4 . 741-03 

-5.041-03  -5.141-03 

-4 

261-04 

52 

0.53 

2.0 

42.0 

41.0 

4 . 43E-43 

-5.391-03  -4.431-03 

*4 

921-04 

52 

0.53 

2.5 

44.3 

45.2 

7. 971-03 

-5.711-03  -4.331-03 

-5 

601-04 

51 

0.52 

3.0 

46.4 

49.3 

7.39B-03 

-6.001-03  -3.671-03 

261-04 

51 

0.51 

3.5 

44  .  1 

93.5 

6 . 671-03 

-6.264-03  -2.441-03 

-4 

494-04 

50 

0.51 

4.0 

49.5 

97.6 

S.42E-03 

-6.471-03  -1.944-03 

-7 

474-04 

51 

0.51 

4.5 

50.6 

101. 

4.43E-03 

-6.644-03  -1.031-03 

-7 

974-04 

51 

0.52 

5.0 

51.2 

106  . 

3.64E-03 

-6.754-03  -5.654-05 

-1 

53 

0.54 

5.5 

51.4 

110. 

2.34S-03 

-6.414-03  9.144-04 

-4 

614-04 

54 

0.55 

4.0 

51.1 

115. 

9.07E-04 

-6.791-03  1.441-03 

-4 

704-04 

9  •J 

0.51 

4.5 

50.2 

119. 

-7.49E-04 

-6.704-03  2.694-03 

-4 

694-04 

60 

0  61 

7.0 

44.7 

124  . 

-2. 604-03 

-6.524-03  3.424-03 

-• 

244-04 

63 

0.64 

7.5 

46.4 

130  . 

-4.67E-03 

-6.244-03  4.004-03 

-7 

734-04 

67 

0.64 

4.0 

43.4 

135. 

-6.96E-03 

-5.464-03  4.434-03 

-4 

934-04 

69 

0.71 

4.5 

39.6 

141  . 

-9. 454-03 

-5.374-03  4.914-03 

-3 

70 

0.71 

9.0 

34.4 

147. 

-1.21E-02 

-4.744-03  5.644-03 

-4 

694-04 

67 

0.64 

9.5 

29 . 1 

152  . 

-1 .4*4-02 

-3.944-03  6.444-03 

-3 

364-04 

61 

0.63 

0.0 

22.3 

157  . 

-1 . 774-02 

-3.064-03  4.444-03 

-2 

034-04 

0.54 
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(•) 

La  t 

Ua 

vo 

Py 

VCor 

Vdrag 

Gaoa 

HA* 

dVo/ 

B 

70.  S 

14.3 

161. 

-2.06E-02 

-1.988-03 

1.178-02 

-6.658-05 

0.42 

0.43 

6 

71.0 

5.13 

164. 

-2.36E-02 

-7.198-04 

1.598-02 

-1.138-05 

0.27 

0.26 

4 

71.5 

-5.51 

165. 

-2.66E-02 

7.498-04 

2.168-02 

-1.34E-OS 

0.10 

0.10 

2 

72.0 

-17.6 

165. 

-2.96E-02 

2.448-03 

2.908-02 

-1.428-04 

-o.io 

-0.10 

0 

72.5 

-31.5 

164. 

-3.2SE-02 

4.378-03 

3.658-02 

-4.668-04 

-0 .34 

-0.35 

B 

73. 0 

-47.2 

160. 

-3.538-02 

6.578-03 

5.038-02 

-1.068-03 

-0.50 

-0.60 

6 

73.5 

-64.4 

154. 

-3.79E-02 

9.008-03 

6.378-02 

-2.078-03 

-0.01 

-0.83 

4 

74.0 

-83.0 

146. 

-4.048-02 

1.168-02 

7.628-02 

-3.568-03 

-1.0 

-1.0 

2 

74.5 

-102. 

137. 

-4.24E-02 

1.448-02 

9.318-02 

-5.628-03 

-1.1 

-1.2 

0 

75.0 

-123. 

127. 

-4.41E-02 

1.738-02 

0.10 

-6.378-03 

-1 . 3 

-1.3 

B 

75.5 

-144. 

116. 

-4.538-02 

2.038-02 

0.12 

-1.168-02 

—1  •  4 

-1.4 

6 

76.0 

-165. 

103. 

-4.588-02 

2.338-02 

0.13 

-1.628-02 

-1.3 

-1.5 

4 

76.5 

-186. 

90.8 

-4.578-02 

2.648-02 

0.14 

-2.158-02 

-1 .  0 

-1.6 

2 

77.0 

-208. 

77.0 

-4 . 478-02 

2.958-02 

0.14 

-2.788-02 

-1.7 

-1.7 

0 

77.5 

-229. 

62.7 

-4.298-02 

3.258-02 

0.15 

-3.518-02 

—  1.7 

-1.6 

B 

70 .0 

-249. 

48.1 

-4.008-02 

3.558-02 

0.15 

-4.358-02 

—  1.0 

-1.6 

6 

7S. 5 

-269. 

33.1 

-3.628-02 

3.848-02 

0.14 

-5.308-02 

—  1.0 

-1.6 

4 

79.0 

-288. 

17.9 

-3.188-02 

4.118-02 

0.13 

-6.358-02 

-l.f 

-1.9 

2 

79.5 

-306. 

2.54 

-2.698-02 

4.388-02 

0.12 

-7.518-02 

-1.0 

-1.9 

0 

80.0 

-323. 

-13.0 

-2.168-02 

4.628-02 

0.11 

-6.778-02 

-l.f 

-1.9 

S 

B0. 5 

-338  . 

-28.7 

-1.628-02 

4.848-02 

0.10 

-0.10 

-1 . 0 

-1.9 

6 

81.0 

-351 . 

-44.5 

-1.088-02 

5.048-02 

9.068-02 

-0.11 

—l.f 

-1.9 

4 

81 .5 

-362. 

-60.3 

-5.528-03 

5.218-02 

7.638-02 

-0.13 

—1  •  9 

-1.9 

2 

82 .0 

-371 . 

-76.1 

-5.818-04 

5.358-02 

6.728-02 

-0.14 

—  1.0 

-1.9 

0 

82.5 

-378. 

-91.9 

3.861-03 

5. 458-02 

5.798-02 

-0.16 

—1.0 

-1.9 

8 

83.0 

-382  . 

-107. 

7.678-03 

5.518-02 

5.108-02 

-0.17 

—1 . 0 

-1.9 

6 

83.5 

-383  . 

-123. 

1.088-02 

5. 538-02 

4.658-02 

-0.19 

—1 . 0 

-1.9 

4 

84.0 

-382. 

-138. 

1.358-02 

5.528-02 

4.408-02 

-0.20 

—1  •  0 

-1.9 

2 

84.5 

-378. 

-154. 

1.578-02 

5.478-02 

4.338-02 

-0.22 

—  1.0 

-1.9 

0 

85. 0 

-373. 

-169. 

1.758-02 

5.398-02 

4.408-02 

-0.23 

-1.0 

-1.9 

s 

85.5 

-365. 

-185. 

1.908-02 

5.288-02 

4.628-02 

-0.25 

-1.0 

-1.9 

6 

86.0 

-354  . 

-200. 

2.018-02 

5.148-02 

4.948-02 

-0.26 

-1.0 

-1.9 

4 

86.5 

-342. 

-216. 

2.108-02 

4.968-02 

5.358-02 

-0.28 

-1.0 

-1.9 

2 

87.0 

-328  . 

-231. 

2.168-02 

4.768-02 

5.648-02 

-0.30 

-2.0 

-1.9 

0 

87.5 

-311. 

-247. 

2.218-02 

4.538-02 

6.368-02 

-0.33 

-2.0 

-1.9 

Table  C.l  Constant  Longitude  Simulation  Data  -  corn. 


DATA  fUOH  SECOBD  SIMULATED  SATELLITE  EUR  ( THESE  OEBITS) 


Tiaa 

Lat 

Lon 

Un 

Vn 

Vdrif 

VCoc 

FT 

Onoa 

MAv 

.4479 

30 . 00 

23.24 

0.00 

0.00 

0.0008*00 

0.0008*00 

0 . 0008*00 

0408+00 

0 . 0008  +  00 

.4501 

30 . 53 

23.25 

0.00 

0.00 

0.0008*00 

0.0008*00 

0.0008*00 

4498+00 

0.0008+00 

.4524 

31.07 

23.21 

0.00 

0 .00 

0 .0008*00 

0 . 0008*00 

0.0008*00 

0408+00 

0 . 0008  +  00 

.  4546 

31 .60 

23.14 

0 .00 

0 .00 

0 . 0008*00 

0.0008*00 

0 . 0008*00 

0448+00 

0.0008*00 

.  4564 

32.13 

23.15 

0 .00 

0.00 

0 .0008*00 

0.0008*00 

0.0008*00 

4008+00 

0.0008*00 

.4590 

32.67 

23  11 

0 .00 

0 .00 

0 .0008*00 

0.0008*00 

0.0008*00 

0408+00 

0.0008*00 

.4612 

33.20 

23.04 

0 .00 

0.00 

0.0008*00 

0.0008*00 

0.0008*00 

•408+00 

0 . 0008*00 

.4635 

33.73 

23.05 

0 . 00 

0.00 

0.0008*00 

0 . 0008*00 

0.0008*00 

•048+00 

0 . 0008*00 

.4657 

34.27 

23.01 

0 . 00 

0.00 

0  .  0008*00 

0.0008*00 

0.0008*00 

•408+00 

0.0008*00 

.4479 

34.80 

22.94 

0.00 

0 .00 

0 .0008*00 

0.0008*00 

0.0008*00 

•••8*00 

0 . 0008*00 

.4701 

35.33 

22.95 

0.00 

0 .00 

0 . 0008*00 

0.0008*00 

0.0008*00 

•048*00 

0.0008*00 

.4724 

35.87 

22.91 

0.00 

0.00 

0.0008*00 

0 . 0008*00 

0.0008*00 

4008*00 

0.0008*00 

.  4740 

36 . 40 

22.41 

0.00 

0.00 

0 . 0008*00 

0.0008*00 

0.0008*00 

•408*00 

0.0008*00 

.476* 

36.93 

22.15 

0.00 

0 .00 

0.0008*00 

0.0008*00 

0.0008*00 

0408*00 

0.0008*00 

.4790 

37.47 

22.41 

0 . 00 

0.00 

0 . 0008*00 

0.0008*00 

0.0008*00 

0908*00 

0.0008*00 

.4012 

38.00 

22.74 

-3.43 

35.34 

-0.4308-02 

0.2618-03 

0.9308-02 

-0 

1368-05 

0.0008*00 

.  4035 

38.53 

22.75 

-4.00 

35.32 

-0 .4328-02 

0.3168-03 

0.9338-02 

-0 

1498-05 

0.1428-04 

.4*57 

39.07 

22.71 

-4.54 

35.29 

-0 .4338-02 

0.3748-03 

0.9368-02 

-4 

2528-05 

0.2268-04 

.4*79 

39.60 

22.64 

-5.17 

35.24 

-0.4358-02 

0.4338-03 

0.9348-02 

-4 

3248-05 

0.3168-04 

.4901 

40.13 

22.65 

-5.77 

35.14 

-0.8358-02 

0.4958-03 

0.9408-02 

-4 

4168-05 

0.4128-04 

.4924 

40.67 

22.61 

-6 . 39 

35.09 

-0 .8368-02 

0.5608-03 

0.9418-02 

-4 

5208-45 

0.5168-04 

.4946 

41  .  20 

22.54 

-7.03 

34.99 

-0.8358-02 

0.6248-03 

0.9428-02 

-4 

6418-05 

0.6278-04 

.496* 

41.73 

22.55 

-7.69 

34.47 

-0.4358-02 

0.7008-03 

0.9428-02 

-4 

7018-05 

0.7448-04 

.4990 

42.27 

22.51 

-4.37 

34.72 

-0.4358-02 

0.7758-03 

0.9418-02 

-4 

9438-05 

0.4798-04 

.  5012 

42.80 

22.44 

-9.07 

34.55 

-0.8348-02 

0.4548-03 

0.9398-02 

-4 

1138-04 

0.1018-03 

.5035 

43.33 

22.45 

-9.40 

34.36 

-0 .4338-02 

0.9378-03 

0.9368-02 

-4 

1348-04 

0.1128-03 

.  5057 

43.87 

22.41 

-10 .55 

34.14 

-0.8318-02 

0.1028-02 

0.9328-02 

-4 

1598-44 

0.1208-03 

.  5079 

44.40 

22.34 

-11.31 

33.91 

-0.8298-02 

0.1118-02 

0.9248-02 

-4 

1468-44 

0.1268-03 

.5101 

44.93 

22.35 

-12.06 

33.68 

-0.8268-02 

0.1208-02 

0.9248-02 

-4 

2158-44 

0.1298-03 

.5124 

45.47 

22.31 

-12.42 

33.43 

-0.8238-02 

0 .1298-02 

0.9208-02 

-4 

2478-04 

0.1298-03 

.5146 

46.00 

22.24 

-13.55 

33.19 

-0.8198-02 

0.1348-02 

0.9168-02 

-4 

2428-44 

0.1278-03 

.5168 

46 . 53 

22.25 

-14.27 

32.95 

-0.8138-02 

0.1478-02 

0.9138-02 

-4 

3198-44 

0.1238-03 

.5190 

47.07 

22.21 

-14.96 

32.72 

-0.8078-02 

0.1558-02 

0.9108-02 

-4 

0.1178-03 

.5212 

47.60 

22.14 

-15.61 

32.50 

-0.7998-02 

0.1648-02 

0.9098-02 

-4 

3968-04 

0.1098-03 

.5235 

48.13 

22.15 

-16.23 

32.30 

-0.7908-02 

0.1728-02 

0.9048-02 

-4 

4368-44 

0.1028-03 

.5257 

48.67 

22.11 

-16 . 79 

32.11 

-0.7798-02 

0.1798-02 

0.9048-02 

-4 

4758-04 

0.9578-04 

.5279 

49.20 

22.04 

-17.32 

31.92 

-0.7668-02 

0.1478-02 

0.9078-02 

-4 

5158-44 

0.9168-04 

.5301 

49.73 

22.05 

-17.41 

31.75 

-0.7528-02 

0.1948-02 

0.9068-02 

-4 

0.4918-04 

.5324 

50.27 

22.01 

-14.26 

31.57 

-0.7378-02 

0.2008-02 

0.9048-02 

-4 

5948-04 

0.4438-04 

.5346 

50  .SO 

21.94 

-14.66 

31.40 

-0.7208-02 

0.2078-02 

0.9008-02 

-4 

6338-44 

0.4908-04 

.5368 

51.33 

21.95 

-19.02 

31.21 

-0 . 7018-02 

0.2128-02 

0.4938-02 

-4 

6708—04 

0.9148-04 

.5390 

51.87 

21.91 

-19.33 

31.03 

-0.6018-02 

0.2148-02 

0.4448-02 

-4 

7468-04 

0.9538-04 

.5412 

52 . 40 

21.44 

-19.60 

30.83 

-0.6588-02 

0.2238-02 

0 .4708-02 

-0 

7408-04 

0.1008-03 

.5435 

52 .93 

21  .  45 

-19.42 

30.62 

-0.6348-02 

0.2248-02 

0.4538-02 

-4 

7718-04 

0.9948-04 

.5457 

53.47 

21.41 

-19.96 

30.42 

-0.6098-02 

0.2328-02 

0 .4328-02 

-4 

7978-04 

0.4748-04 

.  5479 

54  .00 

21.74 

-19.99 

30.26 

-0 . 5838-02 

0.2348-02 

0.4108-02 

-4 

4168-44 

0.6078-04 

.5501 

54.53 

21  .  75 

-19.49 

30.17 

-0 .5588-02 

0.2358-02 

0.7448-02 

-4 

4248-94 

0.1978-04 

.5524 

55.07 

21.71 

-19.53 

30.18 

-0.5348-02 

0.2348-02 

0 . 7678-02 

-4 

4148-04 

-0.3638-44 

.5546 

55-60 

21.64 

-19.16 

30.32 

-0.5138-02 

0.2318-02 

0 . 7508-02 

-4 

7958-44 

-4.1098-03 

.5568 

56.13 

21.65 

-14.46 

30.63 

-0.4958-02 

0.2248-02 

0.7348-02 

-4 

7538-44 

-0.2008-03 

.5590 

56.67 

21.61 

-17.44 

31.14 

-0 . 4428-02 

0.2138-02 

0.7328-02 

-4 

6498-44 

-0.3138-03 

.5612 

57 . 20 

21.54 

-16.14 

31.49 

-0 .4758-02 

0.1998-02 

0.7358-02 

-4 

6438-44 

-4.4338-01 

.5635 

57.73 

21  .55 

-14.52 

32.92 

-0 .4748-02 

0.1798-02 

0 . 7448-02 

-4 

4958-44 

-4.6148-03 

.5657 

58.27 

21.51 

-12.44 

34.25 

-0 . 4408-02 

0.1558-02 

0.7728-02 

-0 

3748-04 

-0 . 7978-03 

.  5679 

51.80 

21.44 

-10.09 

35.45 

-0 . 4938-02 

0.1258-02 

0.4068-02 

-0 

2498-04 

-0.9468-03 

.  5701 

59.33 

21.45 

-7.36 

37.70 

-0 .5128-02 

0.9148-03 

0 .4498-02 

-0 

1368-44 

-0.1198-02 

.5724 

59.87 

21.41 

-4.32 

39.40 

-0 . 5378-02 

0 . 5348-03 

0.9018-02 

-0 

4778-05 

-0 .1318-02 

.  9479 

30.00 

0.74 

0.00 

0 .00 

0 . OOOE+OO 

0.0008*00 

0.0008+00 

0 

•008*00 

0 . 0008+00 

.  9501 

30.53 

0.75 

0.00 

0  .00 

0 . 000E+00 

0 . 0008*00 

0 . 0008+00 

0 

0008*00 

0 . 0008+00 

.  9524 

31.07 

0.71 

0.00 

0  .00 

0 . 0008*00 

0 . 0008*00 

0 . 0008+00 

0 

•008+00 

0.0008+00 

.9546 

31.60 

0 . 64 

0.00 

0  .  00 

0  .  0008*00 

0 . 0008*00 

0.0008+00 

0 

0008+00 

0.0008+00 

.9564 

32.13 

0.65 

0.00 

0  .00 

0 .0008*00 

0 . 0008*00 

0.0008+00 

0 

0008+00 

0 . 0008  +  00 

Table  C.2  Constant  Latitude  Simulation  Data 
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DATA  F10M  SIC01D  SIMULATED  SATELLITE  SUE  ( THESE  OIBITS) 


-  continued  fro*  previous  p»g*  - 


Tias 

Lst 

Lon 

On 

Vn  7dr*9 

VCor 

fy  Gsoa 

5.9500 

32.47 

0.61 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9412 

33.20 

0 . 50 

0.00 

0.00  0.0001*00 

0 . 0001*00 

0.0001+00  0.0001+00 

5.9435 

33.73 

0.55 

0.00 

0.00  0.0001*00 

0 . 0001*00 

0.0001+00  0.0001+00 

5.9457 

34.27 

0 .51 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9479 

34.00 

0.40 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9701 

35.33 

0.45 

0.00 

0.00  0.0001*00 

0.0001*00 

0. 0001+00  0.0001+00 

5.9724 

35.07 

0.41 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9744 

34.40 

0.30 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9740 

34.93 

0.35 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9790 

37.47 

0.31 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

5.9412 

30.00 

0.20 

-5.65 

35.01  -0.7741-02 

0.5051-03 

0.0451-02  -0.3701-05 

5.9035 

34.53 

0.25 

-6.54 

35.64  -0.7661-02 

0.5951-03 

0.0431-02  -0.5051-05 

5.9057 

39.07 

0.21 

-7.37 

35.47  -0.7501-02 

0.6001-03 

0.0611-02  -0.6531-05 

5.9479 

39.40 

0.10 

-0.15 

35.29  -0.7501-02 

0.7611-03 

0.0601-02  -0.0151-05 

5.9901 

40.13 

0.15 

-0.09 

35.11  -0.7421-02 

0.0391-03 

0.0501-02  -0.9071-05 

5.9924 

40.47 

0.11 

-9.57 

34.92  -0.7331-02 

0.9131-03 

0.0561-02  -0.1171-04 

5.9944 

41.20 

0.00 

-10.22 

34.73  -0.7241-02 

0.9031-03 

0.0541-02  -0.1361-04 

5.9940 

41.73 

0.05 

-10.02 

34.52  -0.7151-02 

0.1051-02 

0.0531-02  -0.1551-04 

5.9990 

42.27 

0.01 

-11.37 

34.30  -0.7061-02 

0 . 1111-02 

0.0511-02  -0.1741-04 

4.0012 

42.00 

-0.02 

-11.00 

34.07  -0.6971-02 

0.1171-02 

0.0491-02  -0.1941-04 

4.0035 

43.33 

-0.05 

-12.34 

33.02  -0.6001-02 

0.1231-02 

0.0471-02  -0.2131-04 

4.0057 

43.47 

-0.09 

-12.76 

33.57  -0.6701-02 

0.1201-02 

0.1461-02  -0.2321-04 

4.0079 

44.40 

-0.12 

-13.12 

33.31  -0.6601-02 

0.1331-02 

0.0441-02  -0.2501-04 

4.0101 

44.93 

-0.15 

-13.43 

33.06  -0.6501-02 

0 .1301-02 

0.0421-02  -0.2671-04 

4.0124 

45.47 

-0 .19 

-13.60 

32.01  -0.6471-02 

0.1411-02 

0.0401-02  -0.2021-04 

4.0144 

44.00 

-0.22 

-13.06 

32.50  -0.6341-02 

0. 1451-02 

0.0391-02  -0.2951-04 

4.0140 

44.53 

-0.25 

-13.97 

32.36  -0.6211-02 

0.1471-02 

0.0371-02  -0.3051-04 

4.0190 

47.07 

-0 . 29 

-14.01 

32.16  -0.6061-02 

0.1491-02 

0.0351-02  -0.3131-04 

4.0212 

47.60 

-0.32 

-13.96 

31.99  -0.5901-02 

0.1501-02 

0.0331-02  -0.3161-04 

4.0235 

44.13 

-0.35 

-13.61 

31.05  -0.5711-02 

0.1491-02 

0.0311-02  -0.3161-04 

4.0257 

40.47 

-0 . 39 

-13.56 

31.76  -0.5521-02 

0. *401-02 

0.0301-02  -0.3101-04 

4.0279 

49.20 

-0.42 

-13.10 

11.72  -0.5321-02 

0.1451-02 

0.0291-02  -0.2991-04 

4.0301 

49.73 

-0.45 

-12.64 

31.73  -0.5131-02 

0.1411-02 

0.0291-02  -0.2001-04 

4.0324 

50.27 

-0 . 49 

-11.94 

31.02  -0.4941-02 

0.1341-02 

0.0311-02  -0.2541-04 

4.0344 

50.00 

-0.52 

-11.05 

31.90  -0.4771-02 

0.1251-02 

0.0341-02  -0.2221-04 

4.0344 

51.33 

-0.55 

-9.93 

32.24  -0.4621-02 

0.1131-02 

0.0391-02  -0.1031-04 

4.0390 

51.07 

-0 . 59 

-0.57 

32.61  -0.4501-02 

0.9031-03 

0.0471-02  -0.1391-04 

4.0412 

52.40 

-0.62 

-6.93 

33.10  -0.4411-02 

0.0001-03 

0.0501-02  -0.9251-05 

4.0435 

52.93 

-0.65 

-4.99 

33.72  -0.4361-02 

0.5701-03 

0.0721-02  -0.4001-05 

4.0457 

53.47 

-0.69 

-2.76 

34.49  -0.4351-02 

0.3191-03 

O.0O9B-O2  -0.1531-05 

4.0479 

54.00 

-0.72 

-0 .29 

35.42  -0.4301-02 

0.2061-04 

0.9001-02  -0.1711-07 

4.0501 

54.53 

-0.75 

2 . 39 

36.53  -0.4421-02 

-0.2901-03 

0.9201-02  -0.1191-05 

4.0524 

55.07 

-0 . 79 

5.25 

37.04  -0.4401-02 

-0.6321-03 

0.9501-02  -0.5051-05 

4.0544 

55.40 

-0.02 

0 .25 

39.36  -0.4551-02 

-0 . 9961-03 

0.9711-02  -0.1471-04 

4.0544 

54.13 

-0.05 

11.36 

41.13  -4.4621-02 

-0.1301-02 

0.9911-02  -6.2051-04 

4.0590 

54.47 

-0 .09 

14.56 

43.16  -0.4691-02 

-0.1771-02 

0.1011-01  -0.4701-04 

4.0412 

57.20 

-0.92 

17.01 

45.40  -0.4751-02 

-0.2101-02 

0.1031-01  -0.7301-04 

4.0435 

57.73 

-0.95 

21 . 10 

40.13  -4.4701-02 

-0.2591-02 

0.1041-01  -0.1051-03 

4.0457 

50.27 

-0 . 99 

24.40 

51.11  -0.4001-02 

-0.3021-02 

0.1051-01  -0.1431-03 

4.0479 

50.00 

-1 .02 

27.72 

54.40  -0.4001-02 

-4.3451-02 

0.1061-01  -0.1001-03 

4.0701 

59.33 

-1.05 

31.09 

50.00  -4.4701-02 

-0.3091-02 

0.1071-01  -0.2421-03 

4.0724 

59.07 

-1.09 

34.49 

61.00  -0.4751-02 

-0.4331-02 

0.1041-01  -0.3041-03 

7.4479 

30.00 

-21.72 

0.00 

0.40  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4501 

30.53 

-21.75 

0.00 

0.00  0.0001*00 

0 . 0001*00 

0.0001+00  0.0001+00 

7.4524 

31.07 

-21.79 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4544 

31.40 

-21 .02 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4544 

32.13 

-21 .05 

0.00 

0.04  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4590 

32.67 

-21.09 

0.00 

0.04  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4412 

33.20 

-21.92 

0.00 

0.00  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4435 

33.73 

-21.95 

0.00 

0.04  0.0001*00 

0.0001*00 

0.0001+00  0.0001+00 

7.4457 

34.27 

-21.99 

0.00 

0.00  0.0001*00 

0 . 0001*00 

0.0001+00  0.0001+00 

7.4479 

34.00 

-22.02 

0.00 

0.00  0.0001*00 

0 . 0001*00 

0.0001+00  0.0001+00 

Table  C.2  Constant  Latitude  Simulation  Data  -  corn. 


u« 

0.0001+00 
0 . 0401+00 
0  .  0001+00 
0.0001+00 
0 . 0001+00 
0.0001+00 
i.simti 
0.0001+00 
0  .  0001+00 
0.0001+00 
O.tllEitt 
0.II1B-M 
0 . 9771-04 
O.IHE-O 
0.1031-03 
0 . 1071-03 
0 . 1121-03 
0.1171-03 
0.1241-03 
0  .USE-43 
0.1341-03 
0.13CE-03 
0 . 13SE-A3 
• .1321-03 
0.1261-03 
0  .  USE— 43 
0.147E-43 
0 . S39E-44 
0 . 7441-04 
0. SOCK-44 
0 . 3511-04 
0 . 7171—03 
0 . 2571-04 
0. *411-04 
4. 1001-43 
4. 1011-03 
0 . 221E-43 
0.202X-03 
0 .3731-03 
0 . 4471-03 
0 .3701-03 
0.7041-43 
0 . 0541-03 
0.1031-43 
6. 1241-42 
0.1501-42 
0.1001-42 
0.21C1-S2 
4.2551-42 
0.2001-02 
0.3451-02 
0.3431-02 

0 .0001*00 
0 . 0001+00 
0 . 0001+00 
0 . 0001+00 
0 . 0001*00 
0.0001*00 
0.0001*00 
0.0001*00 
0 . 0001*00 
0 .0001*00 
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DATA  FROM  S8COBD  SIMULATED  SATELLITE  RUB  (THREE  ORBITS) 


Tins 

Let 

Lon 

Ua 

7.4701 

35.33 

-22.05 

0.00 

7.4724 

35.87 

-22.09 

0.00 

7.4746 

36.40 

-22.12 

0.00 

7.4768 

36.93 

-22.15 

0.00 

7.4790 

37.47 

-22.19 

0.00 

7.4612 

38.00 

-22.22 

-4.45 

7.4835 

38.53 

-22.25 

-5.15 

7.4857 

39.07 

-22.29 

-5.76 

7.4879 

39.60 

-22.32 

-6.29 

7.4901 

40.13 

-22.35 

-6.77 

7.4924 

40.67 

-22.39 

-7.19 

7.4946 

41.20 

-22.42 

-7.56 

7.4968 

41.73 

-22.45 

-7.89 

7.4990 

42.27 

-22.49 

-8 .18 

7.5012 

42.60 

-22.52 

-6.44 

7.5035 

43.33 

-22.55 

-1.67 

7.5057 

43.87 

-22.59 

-4.83 

7.5079 

44 . 40 

-22.62 

-4.91 

7.5101 

44.93 

-22.65 

-4.89 

7.5124 

45.47 

-22.69 

-8.74 

7.5146 

46.00 

-22.72 

-8.45 

7.5166 

46.53 

-22.75 

-7.97 

7.5190 

47.07 

-22 . 79 

-7.29 

7.5212 

47.60 

-22.82 

-6.37 

7.5235 

48.13 

-22.85 

-5.19 

7.5257 

48.67 

-22.89 

-3.75 

7.5279 

49.20 

-22.92 

-2.06 

7.5301 

49 . 73 

-22.95 

-0.12 

7.5324 

50 . 27 

-22.99 

2.08 

7.5346 

50.80 

-23.02 

4.55 

7.5366 

51 .33 

-23.05 

7.29 

7.5390 

51.87 

-23.09 

10.31 

7.5412 

52.40 

-23.12 

13.64 

7.5435 

52.93 

-23.15 

17.28 

7.5457 

53.47 

-23.19 

21.19 

7.5479 

54 . 00 

-23.22 

25.33 

7.5501 

54.53 

-23  .  25 

29.63 

7.5524 

55.07 

-23 . 29 

34.05 

7.5546 

55.60 

-23.32 

38 . 54 

7.5568 

56.13 

-23 . 35 

43.05 

7.5590 

56.67 

-23.39 

47.54 

7.5612 

57.20 

-23.42 

51 . 94 

7.5635 

57.73 

-23.45 

56  .  22 

7.5657 

58 . 27 

-23.49 

60.36 

7.5679 

56 .80 

-23.52 

64.37 

7.5701 

59.33 

-23.55 

68 .27 

7.5724 

59.87 

-23.59 

72.08 

-  continued  Iron  previous  pig* 

Vn  Vd rag  VCor 

0.00  0.0008+00  0 . OOOE+OO 

0.00  0 . 0008+00  0. OOOE  +  OO 

0.00  0. OOOE+OO  0. 000E+00 

0.00  0. OOOE+OO  0. OOOE+OO 

0.00  0. OOOE+OO  0. 000B+00 

3  2.01  -0 . 7468-02  0.3718-03 

31.00  -0 . 7368-02  0.441B-03 

31.76  -0 . 7278-02  0.504E-03 

31.65  -0.717E-02  0.5608-03 

31.57  -0 . 700B-02  0. 6111-03 

31.52  -0.699B-02  0.6568-03 

31.49  -0.690E-02  0.700E-03 

31.49  -0 . 6018-02  0.740E-03 

31.51  -0 . 672E-02  0.770E-03 

31.55  -0 . 663B-02  0.015E-03 

31.62  -0 . 653E-02  0.049E-03 

31.72  -0 . 643E-02  0. 0708-03 

31.05  -0.6348-02  0.900E-03 

32.01  -0 . 62SE-02  0.9128-03 

32.21  -0.610E-02  0 . 9108-03 

32.46  -0 . 611E-02  0.092E-03 

32.76  -0 . 607E-02  0.054E-03 

33.12  -0 . 604E-02  0.794E-03 

33.54  -0 . 6038-02  0.707E-03 

34.03  -0 . 605E-02  0.591E-03 

34.62  -0 . 609E-02  0.445E-03 

35.33  -0 . 615E-02  0.269B-03 

36.19  -0.624E-02  0.607E-04 

37.22  -0 . 636S-02  -0.100E-03 

30.46  -0.6498-02  -0.4S4E-03 
39.93  -0. 6651-02  -0.764E-03 
41.67  -0 . 6  0  38-0  2  -0.1118-02 

43.73  -0 . 70  4E-0  2  -0.150E-02 

46.12  -0.7208-02  -0.194E-02 
40.04  -0 . 75  2E-0  2  -0.241E-02 
51.02  -0. 7778-02  -0.291B-02 
55.05  -0. 7991-02  -0.34SE-02 

50.46  -0.010E-02  -0.400E-02 
62.03  -0.0328-02  -0.456E-02 

65.73  -0 . 039E-02  -0.514E-02 

69.51  -0. 0378-02  -0. 5718-02 
73.35  -0 . 024E-02  -0.629E-02 

77.22  -0 . 790E-02  -0.605E-02 
01.09  -0.7598-02  -0. 7418-02 
04.99  -0. 7071-02  -0.795B-02 
00.91  -0 . 644E-02  -0. 0498-02 
92.00  -0.5708-02  -0.9028-02 


Py 

Ceoa 

HAv 

0 .0008+00 

0.0008+00 

0.0008+00 

0 . 0008+00 

0.0008+00 

0.0008+00 

0 . 0008  +  00 

0.0008+00 

0.0008+00 

0 . 0008+00 

0.0008+00 

0.0008+00 

0 . 0008+00 

0.0008+00 

0.0008+00 

0.6388-02 

-0.2308-05 

0.0008+00 

0.6358-02 

-0.3138-05 

0.7638-04 

0.8338-02 

-0.3998-05 

0.6238-04 

0.8328-02 

-0.4868-05 

0.4748-04 

0.8338-02 

-0.5738-05 

0.3368-04 

0.8358-02 

-0.6548-05 

0.2058-04 

0.8388-02 

-0.7428-05 

0.8028-05 

0.0438-02  -0.0238-05  -0.4058-05 
0.0408-02  -0.9028-05  -0.1598-04 
0.0558-02  -0.9798-05  -0.2018-04 
0.0648-02  -0.1058-04  -0.4168-04 
0.0748-02  -0.1118-04  -0.5698-04 
0.0058-02  -0.1158-04  -0.7428-04 
0.0978-02  -0.1178-04  -0.9408-04 
0.9118-02  -0.1158-04  -0.1168-03 
0.9278-02  -0.1108-04  -0.1428-03 
0.9448-02  -0.9948-05  -0.1718-03 
0.9638-02  -0.8468-05  -0.2048-03 
0.9638-02  -0.6598-05  -0.2438-03 
0.1018-01  -0.4468-05  -0.2938-03 
0.1038-01  -0.2388-05  -0.3568-03 
0.1058-01  -0.7328-06  -0.4418-03 
0.1088-01  -0.2518-08  -0.5458-03 
0.1118-01  -0.7748-06  -0.6728-03 
0.1148-01  -0.3768-05  -0.6308-03 
0.1178-01  -0.9658-05  -0.1028-02 
0.1208-01  -0.2018-04  -0.1268-02 
0.1248-01  -0.3588-04  -0.1558-02 
0.1278-01  -0.5868-04  -0.1688-02 
0.1318-01  -0.8998-04  -0.2228-02 
0.1348-01  -0.1318-03  -0.2568-02 
0.1388-01  -0.1638-03  -0.2918-02 
0.1428-01  -0.2468-03  -0.3258-02 
0.1458-01  -0.3228-03  -0.3598-02 
0.1488-01  -0.4098-03  -0.3918-02 
0.1518-01  -0.5098-03  -0.4228-02 
0.1548-01  -0.6218-03  -0.4508-02 
0.1568-01  -0.7428-03  -0.4768-02 
0.1568-01  -0.8738-03  -0.5028-02 
0.1598-01  -0.1018-02  -0.5298-02 
0.1618-01  -0.1178-02  -0.5598-02 
0.1618-01  -0.1338-02  -0.5678-02 


Tabic  C.2  Constant  Latitude  Simulation  Data  -  cont 


Appendix  D  DE  Satellite  Data  Summary 


Much  of  the  data  used  to  compute  the  momentum  forcing  along  the  satellite  tracks 

was  not  explicitly  referenced  in  the  text.  This  appendix  includes  that  data  that,  though 
certainly  important  for  computational  importance,  was  not  included  in  the  discussion. 
Rather  than  take  up  space  in  the  main  body,  it  was  more  prudent  to  include  the  raw  data  in  a 
separate  section.  The  raw  DE2  data  for  the  three  orbits  is  contained  in  the  following  tables 

and  figures: 

Tables  D.l(a  thru  c): 

Raw  data  as  retrieved  from  the  data  base 

Figure  D.  1 : 

Zonal  Neutral  Wind  (Un) 

Figure  D.2: 

Meridional  Neutral  Wind  (Vn) 

Figure  D.3: 

Meridional  Ion  Wind  (Vi) 

Figure  D.4: 

Ion  Temperature  (TO 

Figure  D.5: 

Neutral  Temperature  (Tn) 

Figure  D.6: 

Ion  Density  (ne) 

Figure  D.7: 

Neutral  Density  (nn) 

85 


■AM  DATA  ■CTIItTCO  ««»  OU  IATCU1TI  DATA  IAIK :  OKB JT  B025  DATaBJOlS 
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Tabic  D.l(a)  Archived  DE2  satellite  data  from  orbit  8025. 


DATA  llttllVID  raon  M2  (ATILLlTt  DATA  BAJC:  OIBIT  (02(  DAI»I)01S 


*• . . . 


Nf,fcr^r>r*r,*r‘ffc^r“h‘p*r*p^r“r‘p>*r*«p'r'r»r'r‘r'Pfcr^r.r*r“r.r'r»r*r«r*r*r*r‘^f,*r»r*r* 
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Archived  DE2  satellite  data  from  orbit  8026 


DATA  ICTtllVID  r>OH  Dll  IAT1LUTI  DATA  IAJ1 :  OMIT  *027  DAI«I1I1> 
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Tabic  D.  1(c)  Archived  DE2  satellite  data  from  orbit  8027 


E2  Observed  Zonal  Neutral 


(oas/uj)  un 
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Satcliite-mci  surcd  zonal  neutral  wind  component 


DE2  Observed  Meridional  Neutral  Wind  (Vn) 


Satellite-measured  meridional  wind  component 


E2  Observed 
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Satellite-measured  ion  temperature 


DE2  Observed  Neutral  Density 


Figure  D.7  Satellite-measured  neutral  density  (O  +  N2) 
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